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FOREWORD 


A study was conducted for the NASA AMES Research Center by the Vought 
Corporation to develop improved methodologies for predicting the propulsive 
induced aerodynamics of V/STOL aircraft in Transition/STQL flight. The study 
was performed under NASA Contract WAS 2-11156 with Mr. Richard Christiansen of 
NASA AMES as contract monitor. The Vought efforts in this program were 
accomplished under the direction of Mr. T. D. Beatty who was the Principal 
Investigator for this contract. He was assisted on the contract by Mr. m. K. 
Worthey. Both personnel are from the Flight Technologies directorate of the 
Vought Corporation. 

The authors are particularly indebted to Mrs. D. L. Lewis and Mr. J. W. 
McCharen for their support in the programming of the computer code. 

This report consists of two volumes. The technical discussion of the 
methodology, verification of the techniques and conclusions and recommendations 
are presented in Volusc I. Volume 11 is a detailed user's manual for tne 
computer code developed. 
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1.0 INTRODUCTION ANU SUMMARY 


The capability or predicting the propulsive induced forces and moments on 
a V/STQL aircraft in the transition or STOL region of flight is very important 
to the aircraft designer. These effects will strongly influence the location 
of the jet nozzles and possibly the required thrust of the engines. The 
Vought V/STol- aircraft propulsive effects i YAPE ) program described in this 
report is a computerized method for calculating these propulsive induced 
effects. 

Volume I of this report describes tr.e technical approach and Volume II is 
a description of the data input to the computer code, and the code itself. 

The VAPE program is a single integrated modular computer program for the 
prediction of propulsive induced effects. This rather large computer program 
(over 190 subroutines) consists of four basic components: (1) nacelle inlet 
analysis module, (2) three dimensional lifting potential flow module with 
viscous effects, (3) three jet modules, and (4) a nacelle force and moment 
module. The transmission of data from one component to another is automatic 
with little or no handwork involved. However, the anility to save results for 
a multi >step analysis is also provided. 

In normal usage the body geometry is input to the HESS three dimensional 
potential flow program and the nacelle geometry is input to the inlet analysis 
method. The velocities on the inlet face as predicted by the inlet analysis 
code are then transferred to the HESS program as boundary conditions on the 
inlet face. The HESS program then accesses one of the jet methods to ootain 
the influence velocities due to the jets. The HESS program is then run to 
completion to obtain the flow field about the vehicle. 

The basic framework of the VAPE program is given in Figure 1-1. The 
program consists of a small executive program that controls the activities of 
the other program components. Each of the major components of the VAPE 
program is actually a complete program in itself. They itave been combined 
with the use of the executive main program to give one new and rather large 
program. The discussions throughout this manual will frequently use the term 
"program" in a rather loose sense when referring to any one of the major 
components. 

Each of the major program components except the inlet force and moment 
module can be used as though they are stand-alone programs with all data input 
at the time of execution. The programs are also capable of saving data in 
forms required by the various other programs. By use of appropriate data-unit 
definition cards, any of these data sets can be saved between machine runs. 
This increases the flexibility of the system and permits the user, for 
example, to examine the printed output from one program before continuing 
execution. 

The amount of output produced by each program component depends upon the 
options selected by the user. The print control flags in each component 
should be selected carefully to avoid massive amounts of unwanted output. 


1-1 



The amount or macnire time iboth CPU and 10 times) will also vary witn 
tne options selecteo. aiso, eacn computer installation has its own algorithm 
for calculating the CPU and 10 times. The user will, therefore, have to run 
several jobs using the different program options to determine tne typical 
solution times for his particular installation. In general, the times for the 
potential flow ana bounaary layer programs in tne inlet analysis routine will 
be relatively short, and the HESS 3-0 program much longer. 

In summary, a computerized preoiction metnod for propulsive induced 
forces and moments in transition and short takeoff and landing (STULJ flight 
is presented, ihe method developed was basea on the Vought V/STOL aircraft 
propulsive effects computer program iVaPE). 

The VAPE program is capaole of evaluating: 

o Effects of relative wind about an aircraft 

o Effects of propulsive lift jet entrainment and flow blockage 

o Effects of engine inlet suction on the aircraft flow field 

o Viscous effects on lifting surfaces 

o Determination of engine inlet forces ana moments including inlet 
separation prediction capability 

The effects of relative wind about an aircraft with or without jets 
and/or inlet effects is determined by a very general three-dimensional panel 
method. 

The effects of the propulsive lift jets are determined by one of three 
different jet models which have been extensively modified and/or developed at 
Vought. Some of the major modifications made to the jet models at Vought 
include: 

o Intermediate ground effects calculations 

o Calculations of the flow field in the STOL region of flight 

The effects of engine inlet suction on the aircraft is determined by a 
NASA Lewis code for axisynmetric inlets which has been modified and automated 
at Vought. This method will determine the pressures on the inlet face and 
nacelle Inlet lips. The VAPE program will then utilize these pressures to 
determine the ram drag and forces acting on the inlet. Calculations may also 
be done to determine when and where separation occurs on the inlet lip. 

The various options of the VAPE program have been verified by comparisons 
between calculated and experimental values. 

A computer program code was delivered to NASA and made operational on the 
NASA CDC 76QU computer. A technical manual for this program is contained in 
Volume I of this report. 
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2.0 INPUT INSTRUCTIONS 


The VAPE program requires one primary control card followed by the required 
sets of data cards for each program option to De executed. The sets of data 
furnished must be in the order as specified by the option number on the primary 
control card. 


The general scheme used in describing the input data is shown below: 


Colume Code 
Col umn 

Code 

Routine 

Format 


Explanation - 


Routi ne 

Format Explanation 

Column indicates the starting position on the card for 
each data field. 

The “code" gives the FORTRAN name used in the read 
statement by the program. 

"Routine" indicates the subroutine where the data in read. 

The parameter “FORMAT" which is given right under the 
routine name, indicates the FORTRAN format of the data 
read statement field. The parameter 15 would indicate 
that the parameter is an integer in a field that is 5 
columns wide. Integers should be punched on the right 
side of the field (right justified). The parameter FiO.O 
would indicate a floating point number punched with a 
decimal point i.e. , -12.354). The number may be punched 
anywhere in the field indicated irrespective of the 
decimal point location indicated by the format. The 
parameter E12.6 would indicate a floating point number 
punched with a decimal point i.e., 5.0 x 10 b ). On 
input E and F formats are treated identically. 

Free format indicates that no format is used. Data is 
input separated by one or more blanks, a corana, or a 
slash. Embedded blanks are not permitted. 

The description of the input data is given under 
"explanation". 
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2.1 DATA STORAGE 


In order for the VAPE program to function efficiently, a large amount of 
data must oe transferred between subprograms. Some of this aata is only 
temporary and is stored on "scratch" files which are not retained when no 
longer needed. Certain data which is required oy various routines is stored 
on permanent files which can be saved, so that the program can be interrupted 
and then restarted. Table 2.1-1 lists the file names, and a short description 
of the file contents. Table 2.1-2 presents a summary of the files used by 
file number. 
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Table 2.1-1 
Data Storage Files 


ORIGINAL PAGE 18 
OF POOR QUALITY 


FILE 

FILE CREATED 

FILE 

FILE USED 

FILE 

DESCRIPTION 

UNIT 

IN SUBROUTINE 

ALIAS 

IN SUBROUTINE 

ALIAS 


MTI 

£00 


LOCAL* 

IFUTAP 

SCRATCH FILE 


BSETUP 

IRS IDE 

SOLMOR 

DEVTAP 

RHS COLUMN 

MT2 

EDO 


LOCAL* 


SCRATCH FILE 


COMBYN 


LOCAL* 


SCRATCH FILE 


BSETUP 


INP2 


BOUNDARY 
LAYER STRIP 
INFORMATION 


SOLMOR 

NOUT 

LOCAL 


SCRATCH FILE 


HESS 


LOCAL 


SCRATCH FILE 

MT3 

EOD 


LOCAL* 


SCRATCH FILE 


COMBYN 


LOCAL* 


SCRATCH FILE 


VCOM 

NOUT 

AFORM 

NOUT 

ONSET FLOW 
MATRIX 




OLCMFLO 

NOUT 


MT4 

EOO 


LOCAL 


SCRATCH FILE 


LIFT 


VFMLFT 


GEOMETRY 

DATA 


NOLIFT 


VFMNLF 

NTT 



COLSOL 

NW.NSIG 

PKUTTA 

OLCMFLO 

HESS 


SIGMA MATRIX 

MT5 





INPUT 

MT6 





OUTPUT 

MT7 

SOLMOR 

NS I 

LOCAL 


SCRATCH FILE 

MT8 

EOD 


LOCAL* 


SCRATCH FILE 


COMBYN 


LOCAL* 


SCRATCH FILE 


VFMLFT 

N8 

VFMLFT 

N8 

VELOCITY 

MATRIX 


VCOM 

IONSET 

PKUTTA 

IONSET 

KUTTA MATRIX 

MT9 

EOD 




SCRATCH FILE 


COMBYN 




SCRATCH FILE 


AFORM 

NAIJ 

COLSOL 

NIN 

AIU MATRIX 

MT10 

EOD 


LOCAL* 


SCRATCH FILE 


COMBYN 




scratch file 


AFORM 

NRSIDE 

COLSOL 

RHSTAP 

RIGHT HAND 
SIDE MATRIX 

MT11 

EOO 


LOCAL* 


SCRATCH FILE 


VFMNLF 


VCOM 






— 


-J 
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Data Storage Files (Continued) 


FILE 

FILE CHEATED 

FILE 

FILE USED 

FILE 

DESCRIPTION 

UNIT 

IN SUBROUTINE 

ALIAS 

IN SUBROUTINE 

ALIAS 





AFORM 

NT 

SOURCE 




ULCMFLQ 

NT 

VELOCITY 

MATRIX 

MT12 

EOD 


LOCAL* 


SCRATCH FILE 


Kv^Tl ■ 

MT19 

JETOLD 





MT19 



JET DATA 



MT19 




MT13 

EOD 


LOCAL* 


SCRATCH FILE 


VFMLFT 

NO 

VCOM 

ND 

DIPOLES 

VELOCITY 


COLSOL 

NOUT 

LOCAL 


SCRATCH FILE 

MT14 

HINIT 


INPUT 

VCOM 


REGENERATED 




INLENT 

JETOLD 

RECJET 

JET3IN 


INPUT 


COLSOL 

NS1 

LOCAL 


SCRATCH FILE 


SOLMOR 

NS1 

LOCAL 


SCRATCH FILE 

MT15 

COLSOL 

NLTAPE 

SOLMOR 

NLTAPE 

L-MATRIX 

MT16 

COLSOL 

NTT APE 

SOLMOR 

NTTAPE 

T-MATRIX 

MTi7 

EOD 


COMBYN 


GEOMETRY 

DATA 


vfmnlt 


PKUTTA 

DVIJ 

K'JTTA RIGHT 
SIDE 


VFMLFT 





MT18 

SOLMOR 

NU 

HESS 

NBLSUL 

ADDITIONAL 
SIGMA FROM 
SQLnR 

MT19 





SAME AS MT12 

MT20 

1 


EOD 


INPUT 




EOD j 


INPUT 


■H 




OUTPUT 

OPTION 

MT21 

VCOM 




OUTPUT 

OPTION 


JETOLD 

JET3IN 
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Data Storage Files (Continued) 


Flit 

JUT 

FILE CREATED 
IN SU6Ru‘iT IiME 

FILE 

ALIAS 

FILE USED 
IN SUBROUTINE 

FILE 

ALIAS 

DESCRIPTION 

i*i 1 Li 

CUrlriYN 

CGMFLO 


LOCAL* 

LOCAL 


SCRATCH FILE 
SCRATCH FILE 

MT23 

COrt&Yh 


LOCAL* 


SCRATCH FILE 

HT25 

CCM8YN 


TRANSG 


TRArtSrlW ItO 
COORDINATE DATA 
FROM COME 

MT26 

TRANSG 


YISCUS 


TRANSLATED, 

SCALED 

COORDINATE DATA 
FROM TRAN 

MT27 

VISCUS 


GEOMOO 


DT DATA FROM 
VISCUS 

MT28 

sciRa 


GEOMOO 


X, Y, ALPHA DATA 
FROM SCIRCL 

MT29 

GEOMOO 


E00 


BUILT FOR EOD 
PASS 2 

HT34 

CQMBYN 


INLINT 



MT37 





MASS STORAGE 

MT39 

KINIT 


8SETUP 


BOUNDARY LAYER 
INPUT 

MT40 





MASS STORAGE 






NOTE: 

*-Nay be located 
in more than 
one routine but 
member of major 
group. 
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2.2 EXECUTIVE SYSTEM CONTROL CARDS 


This routine controls the order in which the routines within VAPE are executed 
and initializes some parameters for use in the Stockman Program. 

E-l PRIMARY CONTROL CARO (12, 13, IS, IIP, FlC.O, IIP, F10.2, 515, F5.0) 


COLUMN CODE 

1 MTCK 

3 IPROG 


6 IPRG2 


11 MXPTS 


ROUTINE 

FORMAT EXPLANATION 

INITMT Tape assignment flag 

(12) « 0 Normal output 

s 1 Output assignment of scratch files 


INITMT 

(13) 


INITMT 

15 


INITMT 

110 


Selects stage at which computation 
begins, default value is 1 . 

Execution begins with: 

s 1 Geometry module (Stockman) 

= 2 Potential flow routine (Stockman) 

= 3 Combination routine (Stockman) 

= 4 Translator routine (Stockman) 

= 5 Viscous routine (Stockman) 

= 6 GEOMCD routine (Stockman) 

3 7 HESS routine (3-0 HESS) 

Selects stage at which computation 
stops default value is 7 

Execution stops at: 

= 1 Geometry routine (Stockman) 

= 2 Potential flow routine (Stockman) 

3 3 Combination routine (Stockman) 

s 4 Translator routine (Stockman) 

3 5 VISCOUS routine (Stockman) 

3 6 GEOMOD routine (Stockman) 

s 7 HESS routine (30 HESS) 

Maximum number of points that may be 
generated by the geometry routine, 
default value is MXPTS - 180. Large 
numbers of points* requime/excessive 
computation time, thus this limit can be 
used to stop execution if it is not 
desired to use large amounts of 
computation time. 


* Problems where SCIRCL generates 250 points or less have little 

impact on computer run time. Problems where SCIRCL generates 300 or 
more points will effect computer run time and should be avoided if 
possible. 
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21 

TINX 

initmt 

F1C.0 

Maximum amount of computation time that 
may be used by the potential flow 
solution routine, default value is TIMX = 
200 (seconds). Errors in the geometry 
iryut can cause excessive amounts of time 
to be used In the potential flow solution 
routine. Program execution is terminated 
If this limit is exceeded to avoid 
wasting compuia.ion time. 

31 

I PASS 

INITMT 

no 

If I PASS * 1 o IPRG2 * 1 program 
execution stops after potential flow 
solution (Stockman). 




If IPASS * 1 a IPRG2 =/ 2 or 3 program 
skips translator, VISCOUS a GEOMOO 
routines and proceeds directly to 3-D 
HESS program. 

41 

XRI2 

INITMT 

Surface distance increment added to 
stagnation point or inlet highlite to 
start boundary l«yer solution. 
(Recooraend 0.05) 



F10.2 


51 

NT(l) 

INITMT 

15 

Number of on-body points for closed 
end solution. 

56 

NT(2) 

INITMT 

IS 

Number of on-body points for open 
end solution. 

61 

NP 

INITMT 

15 

Total number of cff-body points 

66 

NHUBMX 

INITMT 

15 

Number of last points on hub 

71 

IPLGT 

INITMT 

15 

= 0 No plots 

= 1 Plots (not checked out) 

76 

XRPLOT 

INITMT 

F5.0 

Last X where data is 

Plotted (Recommend XRPLOT * XRI) 


(1) Note that NT(1), NT(2), NP, and NHUBMX should be input only if IPR0G>1. 
These values are output of SCIRCL. 

(2) The parameters from MXPTS through XRPLOT are control values for the 
nacelle inlet analysis module and are not required if IPROG = 7. 
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2.3 STOCKMAN INLET ANALYSIS PROGRAM* 

G£CHET:« .-SOULE 


IA-1 TITLE CARD («A 10) 


COLUMN 

COPE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

Title 

SCIRUL 

0A1U 

Title card 

If Title (1) = Stop, Control transfers 
to main routine, input to SCIRCL is 
complete. 

In-2 CONTROL CARD 

FOR GEOMETRY MODULE (ZA10.4I5) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

IDENT 

SCIRCL 

A10 

Case identifier 

11 

PROG 

SCIRCL 

A10 

Potential flow flag 

£00 Ax i symmetric program 
Used for identification only. 

21 

N06 

SCIRCL 

15 

Flag for potential flow solution 

» 0 All solutions 

s 1 Basic geometry data only 

* Cards 

for this module required only if IPROG = 1 or. cnru C L 
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IA-3 GEOMETRY DISTRIBUTE CCNTRGL CARD (8F10.2) 


COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

ANBDYS 

SCIRCL 

F10.2 

= 1 
= 2 
= 3 

Number of input bodies 
Shroud (cowl) only 
Hub and shroud 

Hub, flow splitter plates and shroud 

11 

CELS 

SCIRCL 

F10.2 


Spacing between points in region of 
interest, DELS affects total number of points, 
solution accuracy, and computer run time. 

21 

DELSMX 

SCIRCL 

F10.2 


Maximum spacing used far from region of 
interest. DELSMX must be less than inlet duct 
flow width 

31 

XRI 

SCIRCL 
F10 .2 


Axial location of point of interest. 

Point at which surface distance equals zero. 


Spacing between points is kept at 
approximately DELS up to XRI , then spacing 
increases to DELSMX far from region of 
interest. XRI should be located in constant 
area section downstream of the assumed fan 
face of the inlet. 

41 ANNSD SCIRCL Number of noise suppression devices or 

F10.2 NSD splitter (can be zero) 
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I A -4 CONTROL CARD FOR FAKES (2015) 


COLUMN 

COPE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NRAKES 

SCIRCL 

15 

Number of axial locations at which data 
across the passage is desired. At least one 
is required @ XRI. 

1 $ NRAKES 5 25 

IA-5 

RAKE DEFINITION CARD (110, 

3F10.5) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NY 

SCIRCL 

110 

Number of points in rake at XRAK. Total 
number of points for all rakes limited at 200 
NY S 200, Y = (YHI-YL0)/NY 

11 

XRAK 

SCIRCL 

F10.5 

Axial location of rake XRAK,-+i s 
XRAKj 

21 

YLO/ 

SCIRCL 

F1C.5 

Y value of first point on the rake 
(must be an off body point) 

31 

YHI 

SCIRCL 

F10.5 

Y value of last point on the rake 
(must be an off body point) 


Note: One rake definition card for each rake: thus NRAKES number of 
cards input here 
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IA-6 BODY CONTROL CARO (6F1U.0) ONE CARO FOR EACH BUOY 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

TYPBDY 

SCIRCL 

F10.0 

Body identification flag See table below 

11 

ANSEG 

SCIRCL 

F10.0 

Number of segments for the particular 
body, see TYPBDY. If TYPBDY = 0, it is the Y 
centerline used for mirroring. 



TABLE I 

Use of Parameter TYPBDY 


Value of 


VALUE of TPYBDY 


AN BOYS 

0 

1 2 3 


1 


shroud 


2 

centerl i ne 

hub shroud 


3 

centerl i ne 

hub *flow splitter shroud 


Analytical Functions 


The inlet surfaces can be described by combining parts of six available 
analytical functions, as shown in Figure 2-1. The coefficients for each analytical 
function used are determined by the program from the coordinate points input. The 
number of coordinate sets needed varies with each analytical function as listed in 
Table II, and shown in Figure 2-1. 

IA-7 SURFACE CURVE IDENTIFICATION CARD al (10X, 7F10.2) 

ROUTINE 

COLUMN CODE FORMAT EXPLANATION 

11 ENREED SCIRCL code indicating type of curve to be 

F10.2 fitted through the given points, it is also 

the exponent for a superellipse when the 
exponent is to be specified. (See Table II 
Below) 

ENREED * 0.0, SCIRCL program will calculate an 
exponent. Read In 5 coordinates (XIN and 
YIN), (fig. 2-l(a) ) 

1.0 < ENREED < 10.0, superellipse with 
exponent equal to ENREED. Input coordinates 
1,2,4, 5. Read In 0.0 for point no. 3. (Point 
no. 3 Is not used, but points 4 and 5 must be 
in the proper columns.) (fig. 2-l(b)) 
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OF POOR QUALITY 1C01. < ENREED < 101C. The flag 1CC0. added 

to the code of desired superellipse is used on 
the shroud to give finer spacing at the 
highlight. The superellipse going into the 
highlight and the one of the topside of the 
highlight should have this flag. 

ENREED = 1.0, is a straight line, input 2 
coordinates (XIN (1), YIN(l), XIN(2), 

Y IN ( 2 ) ) . (fig. 2-l(c)) . 

ENREED = 10. C, special straight line used for 
closed bodies (example - airfoils). This 
straight line starts with large spacing 
(CELSNX) and ends with the small spacing 
(DELS). Input 2 coordinates, (fig. 2-1 (d) ) - 

ENREED = -1.0, fits a lemniscate between a 
straight line and a point. Input is 3 
coordinates, (fig. 2-l(e) ) 

ENREED * -2.0, fits an ellipse between two 
straight lines. Input 4 coordinates. 
Restriction, the two straight lines should be 
at a right angle, therefore superellipse 
routine can be used. (fig. 2-1 ( f ) ) 

ENREED = -3.0, fits a cubic between 2 straight 
lines, input 4 coordinates, (fig. 2-l( g) ) 


21 

Xlh(l) 

SCIRCL 

F10.2 

X value of first input point 

31 

YIN(l) 

SCIRCL 

F10.2 

Y value of first input point 

41 

XIN (2) 

SCIRCL 

F10.2 

X value of 2nd input point 

51 

YIN(2) 

SCIRCL 

F10.2 

Y value of 2nd input point 

61 

X IK (3 ) 

SCIRCL 

F10.2 

X value of 3rd Input point 

71 

YIN (3) 

SCIRCL 

F10.2 

Y value of 3rd input point 
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Figure 2-1. Analytical functions available for describing Inlet surfaces 
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IA-£ SURFACE CURVE IDENTIFICATION CARE » 2 (8F10.2) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

Xlh (4 ) 

SCIRCL 

X value of 4th input point 



F10.2 


11 

YIN (4) 

SCIRCL 

F10.2 

Y value of 4th input point 

21 

XIN(5) 

SCIRCL 

F10.2 

X value of 5th input point 

31 

YIN (5) 

SCIRCL 

F10.2 

Y value of 5th input point 
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TABLE II 

Analytical Functions Available 


Curve 

Code 

ENREED 

Number of 
Points Needed 

Portion 

Used (Figure 2-1) 

Cubic 

3.0 

4 

2 to 3 

El 1 i pse 

2.0 

5 

2 to 4, point 3 is a 
dummy pt 

Lemni scate 

1.0 

3 

2 to 3 

Super Ellipse 

0.0 

5 

2 to 4 

Straight Line 

1.0 

2 

1 to 2 

Super Ellipse 

1.0 to 10.0 

4 

2 to 4, pt. 3 optional 


1001. to 1010. 

4 

2 to 4, extra fine grid 

Straight Line 

10.0 

2 

straight line for closed 
body 


There should be ANSEG number of the above two cards input at this point. 
Additionally, there should be the last three cards repeated for each body (ANBDYS) 
at this point. 


IA-9 SPLITTER BODY DEFINITION CARD »1 (8F10.2) (input only if ANNSD >0) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

DELS 

SUPRD 

F10.2 

Spacing between points on the NSD 
splitter, can be different than DELS on the 
body 

IA-10 

SPLITTER BODY DEFINITION 

CARD #2 (1615) (input only if ANNSD >0) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NSEG 

SUPRD 

15 

Number of segments on the NSD splitter. 

For a thin NSD splitter, segments on the 
topside should begin and end at the same X 
values as corresponding segments or the 
underside. 

6 

NSHJGH 

SUPRD 

15 

Number of the segment on the underside of 
the NSD going into the highlight, equal to 1/2 
of NSEG for thin splitter. Set NSHIGH = NSEG 
for all other splitters. 
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IA-11 SPLITTER SURFnCE CURVE IDENTIFICATION CARDS 110X, 7F10.2) (Input only 
it ANNSL 1 >U) 


The information on the geometry identification ird and tne surface curve 
identification card must be input for the splitter using same definitions and 
card form as given above for these two cards so input ENREEU, XIMli), YlNil) 
through XIU(5), YIN(5) at this point for splitter 

REMARKS, Total number of points for all the bodies should not exceed 400. 

The first straight line on the shroud must be equal in axial lengtn 
to the last straight segment on the nub. If the first straight 
line on the shroud is longer, use two segments. 

Total number of off-boay points must not exceed 200. If there are 
acoustic splitters in the area of a rake, specify ANNSD + 1 rakes 
for axial location of raxe. 
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AXISYHJETRIC POTENTIAL FLOW PROGRAM 

The Axisymmetric Potential Flow Program is the second part of the 
Stockman inlet analysis procedure. The output from SCIRCL (The Geometry 
Module) is used in this sub-program. The results of this sub-program are used 
by the combination routine to determine the flow field. 

If data is being transferred from SCIRCL and NIN i 0 or NIN i 5 then only 
the following card is input for the potential flow program. 

Input data only if I PROG = 2 on card E-l. 

AP-1 POTENTIAL FLOW CONTROL CARD (215) 




ROUTINE 


COLUMN 

CODE 

FORMAT 

EXPLANATION 

1 

NIN 

EOD 

input tape for EOD, normally generated 



15 

by SCIRCL 




= 5 data input on cards (card must follow) 




* 20 data being transfered from SCIRCL 

6 

IEODU 

EOD 

i 0 extra printout for debugging of program 



15 

(normally set 20) 


If NIN = 5, then the following cards must be input. If NIN = 20, then skip 
to combination program card description. 

AP-2 TITLE CARD (5A1Q, A7 , 5X, 110) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

KEDR 

PARTI 

5A10 

title of case. 

51 

CASE 

PART 1 
A7 

case identi *«* 

63 

SEQ1 

PART i 

no 

sequence number cf this card in the input 
scheme 


AP-3 FLAG CARO (2811, 12, 32X, IIP) 


Card columns 1-30 when punched with ary non-zero Integer, activate flags 
that Indicate the following: 
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CCLITiN 

cGoe 

ROUTINE 

FORMAT 


EXPLANATION 

1 

MS 

PARTI 

11 


The number of bodies input. Normally 
set equal to 1. 1 s MS * 5 

2 

NNU 

PARTI 

11 


The number of non-uniform onset flows. 
Normally set equal to 0. 

3 

FLGC3 

PARTI 

11 

-C 

Axi symmetric flow flag. 

No axi symmetric stream-flow solution 
calculated. 




-1 

Axisynmetric streamflow solution is 
calcul a ted 





Normally set equal to 1 

4 

FLG04 

PARTI 

11 

=0 

Cross flow flag. 

No cross flow solution is calculated 




=1 

Cross flow solution is calculated 





Normally set equal to 0 

5 

FLG05 

PARTI 

11 

=0 

Off-body point flag 
No off body points input 




=1 

Off body points are input 





This flag allows the velocity at points 
off the body surface to be determined. 

6 

FLG06 

PARTI 

11 

=0 

Basic data formation flag 
A full case will be done 




=1 

The basic data, i.e. , midpoints, normals 
etc. will be formed and printed. No 
velocities will be calculated. 

7 

aG07 

PARTI 

11 

=0 

Ellipse generator option 
Body coordinates will be input 




=1 

An ellipse is generated using data input 
later. No body coordinates are input 

8 

FLGQ8 

PARTI 

11 


Matrix print flag 


*0 Coefficient matrices are not printed. 

=1 Coefficient matrices will be printed. 

Normally set equal to 0 
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9 FLG09 PARTI 
II 


10 FLG10 PARTI 
II 


11 FLG11 PARTI 
11 


12 FLG12 PARTI 
II 


13 FLG13 PARTI 
II 


14 FLG14 PARTI 
II 


15 FLG15 PARTI 
II 


16 FLG16 PARTI 
II 


Matrix solution flag#l 
=0 Matrix solution set later 

-1 Old SIEOEL method used 

Matrix solution flag #2 

=0 Matrix solution set elsewhere 

=1 Modified SIEOEL method used 

Perturbation velocity flag 
=0 Normal case 

si No onset flow used. Only perturbation 
velocities are calculated. 

Potential matrix solution 

sO Normal case 

si A potential matrix is -solved 

Matrix solution flag #3 

s' 1 Matrix solution set elsewhere 

si Matrix solution Dy triangularization 

Normally set equal to 1. 

Prescribed tangential velocity flag 

sO Normal case 

=1 Tangential velocities are specified 

Strip ring vorticity flag 
sO Normal case 

=1 A vorticity distribution is formulated. 

Axisynmetric uniform flow flag 
sO Normal case 

si Axi symmetric uniform flow solution is 

omi tted 

Normally set equal to 0. 
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17 

FLG17 

PARTI 


Crossflow uniform flow flag 



11 

=0 

Normal case 




=1 

Crossflow uniform flow solution is 
omitted. 





Since FLGC4 is normally = 0 then so is 
FLG17 normally set equal to 0. 

18 

FLG18 

PARTI 


Surface vorticity flag 



11 

=0 

Normal case 




=1 

Surface vorticity is generated. 

19 

FLG19 

PARTI 


Prescribed vorticity Flag 



11 

=0 

Normal case 




=1 

A prescribed vorticity is input 

20 

FLG20 

PARTI 


Total vorticity flag 



11 

=0 

Normal case 




=1 

Total vorticity calculated 

21 

FLG21 

PARTI 


Extra crossflow flag 



11 

=0 

Normal case 




=1 

Extra crossflow option used 

22 

FLG22 

PARTI 


Generated boundary condition flag 



11 

=0 

Normal case 




=1 

Boundary conditions generated 

23 

FLG23 

PARTI 

11 


Ring wing option flag 



*0 

Normal case 




*1 

Ring wing option used 

24 

FLG24 

PAP.T1 


Not used by this program 



11 



25 

FLG25 

PARTI 


Not used by this program 



11 


26 

FLG26 

PARTI 


Not used by this program 



11 
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27 

FLG27 

PARTI 

11 

Not used by this program 

28 

IPUVEL 

PARTI 

11 

0* 

Punch flag 
No data punched 



=1 

Punched output generated 

29 

NOUN 

PARTI 

11 

Not used by this program 

63 

SEQ2 

PARTI 

11 

Sequence number of this case in the input 
stream 

AP-4 

CHORD CARD 

(3F10.0, 32X, I1C) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

CHORD 

PARTI 

F10.0 

Reference chord length used to non- 
dimensional ize x and y coordinates 

11 

IM 

PARTI 

F10.0 

Mach number (Ml < 1.0) use to approximate 
effect of compressibility (Gothert 1 rule) 

21 

TCNST 

PARTI 

F10.0 

This is a constant which is used for the 
value of the tangential velocity if this 
option is desired. 

62 

SEQ2 

PARTI 

F10.0 

Sequence number of card in input stream 

AP-5 

BOOT TRANSFORMATION CARD (5X, 

15. 5F10.0. 16X. 14) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

5 

NN 

BASIC1 

15 

The number of input points on this body. 
NN s 500 

11 

MX 

BASIC1 

F10.0 

A factor used to multiply all 
x-coordi nates. MX is assumed equal to 1 
If no valve is input. 

21 

MY 

BASIC1 

F10.0 

A factor used to multiply all y- 
coordi nates. MY is assumed equal to 1 if 
no value is input. 

31 

THETA 

BASIC1 

F10.C 

An angle (in degrees) through which all 
points of a body are to be rotated about 


the origin In the clockwise direction. 
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41 

ADOX 

BASIC1 

F10.C 


A constant to be added to all x- 
c cord i nates 

51 

ADDY 

EASIC1 

F1C.C 


A constant to be added to all y- 
cocrdinates 

76 

SEQ2 

BASIC! 

14 


Sequence number of this card in input 
stream 

AP-6 

BODY CONTROL 

CARD (3(5X, 

. 15). 

2F10.0, 26X, 14) 

COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

6 

BON 

BASIC1 

15 


Body sequence number. This program will 
handle up to 5 bodies. 

16 

SUBKS 

BASIC1 

15 

=0 

Subcase Flag. 
Normal case 




=1 

Use unmodified coordinates of the 
previous case. 

26 

NLF 

BASIC1 

15 

=0 

Non-lifting flag 

Body is lifting (this is used in special 
option) 




=1 

Body is non-lifting 
(normal case) 

31 

XE 

BASIC1 

F10.C 


Value of major semi -axis for use by 
ellipse generation option. 

41 

YE 

BASIC1 

F10.0 


Value of minor semi-axis for use by 
ellipse generation option 





Note: if XE s YE a sphere will be 
formed. 

76 

SEQ2 

BASIC1 

F10.0 


Sequence number of this card in the input 
stream 


GEOfCTRY DATA CARDS (6F10.0, 16X, 14) 


The body geometry data cards are Included only if the input parameters 
NIN i 0 or 5 and SUBKS = 0 on the body control card. If NIN = 20 then the 
data is read from unit 20. If NIN * 5 and BDN * 0 then the following cards 
contain the x-y coordinates of off-body points instead of x-y geometry data. 
The number of either geometry data point or off-body points must be equal to 
NN. 
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AP-7 X-CCOF.DINATE CARDS (six values per card) 




ROUTINE 


COLUMN 

CODE 

FORMAT 

EXPLANATION 

1 

TXl(l) 

BASIC1 

6F10.0 

x-coordinates of body input from 
leading to trailing edge 

11 

TX1(2) 

BASIC1 

6F10.0 


21 

TX1(3) 

BASIC1 

6F10.C 


etc. 

76 

SEQ2 

EASIC1 

14 

Sequence number of this card in the 
input stream 


AP-8 Y-COORDINATE CARDS (six values per card) (6F10.0, 16X, 14) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

TY1(1) 

BASIC1 

6F10.C 

y-Coordinates of body which correspond 
to the x-values above, y values must be 
positive. 

11 

TY1(2) 



21 

TY1(3) 



etc. 




76 

SEQ2 

BAS I Cl 
14 

Sequence number of these cards in input 
stream 

NOTE: 

Each body input, including the off body points, requires the body 
transformation card, the body control card, and may also require the 
geometry data cards depending on the input flags. This is the 


stopping place for a normal axisymmetric case. The following cards 
are input only if one of the special options is required. 

AP-9 TANGENTIAL VELOCITY DATA CARDS (six values per card) (6F10.0, 16X, 14) 

These cards are input only if FLG14 * 0 and TCNST = 0.0 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

TG(1) 

BASIC1 

Specified tangential velocities at 

11 

TG(2) 

6F10.0 

element midpoints. 
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21 


etc . 


TC(3) 


76 SEC2 SASIC1 Sequence number of these cards in input 

14 stream 

AP-10 NCU4JNIFCRM FLCW CARDS 


This card is input only of NNU * 0. 


COLUMN 

COCE 

ROUTINE 

FORMAT 

EXPLANATION 

6 

NUM 

BASIC2 

15 

Non-uniform flow identification number. 

16 

MSF 

BASIC2 

15 

If MSF = 1 the flow velocities Nq,T 0 
will be be used for the axi symmetric case 
only. 




If MSF = 1 the flow velocities N 0 ,T 0 will 
be used for the cross flow case only. 




If MSF > 1 the flow velocities will be used 
for both axi symmetric and cross flow cases. 

21 

TYPE 

BASIC2 

F10.0 

Flag which specifies the type of input 
flow velocities at each mid-point. If TYPE > 
0.0, the velocities are input as x & y 
components. 




If TYPE = 0.0 the velocities are input as 
normal & tangential components. 




If TYPE <0.0 the automatic generation of the 
flow due to a rotating body is used. 

31 

FG 

BASIC2 

F10.0 

Constant used by the flow generator. 
Type must be less than 0.0. 

76 

SEQ2 

BASIC2 

14 

Sequence number of this card in the 
input stream 


The following cards 

are input only if NNU 4 0 and TYPE 4 -1.0. 

AP-11 

NORMAL VELOCITY CAROS (six values per card) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NG(1) 

BASIC2 

F10.6 

This is either the x or normal velocity 
component depending on the value of type 
above. These values must be in sequence with 


the coordinate data. If the x component is 
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input it is defined as positive to the 
right. If the normal velocity is input it is 
positive if it is to the interior of the 
body. NN-1 values are input. 


76 

SEQ2 

BASIC1 

14 

Sequence number of this card 
cards in input stream 

AP-12 

TANGENTIAL 

VELOCITY CARDS 

(six values per card) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

TGC1) 

BASIC2 

6F1C.0 

This is either the y or tangential 
velocity component depending on the 
value of type above. These values must 
correspond to the NG values above. If the y 
component is input it is defined as positive 
if it is orientated upwards. If the 
tangential velocity is input it is positive 
if the flow field is to the left of the 
vector representing the tangential velocity. 

76 

SEQ2 

BASIC2 

14 

Sequence number of this card in the 
input stream 


2-26 



INPUT FCR COMBINATION PROGRAM 


C-I TITLE CARD (8A1C) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

TITLE 

READS 

8A10 

Title of this case 

If titled) input as end control returns 
to main program - input for combination 
program complete 


C-2 COMBINATION ROUTINE CONTROL CARD (2015) 


COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

NID 

READS 

15 


Number of ID cards from EOD output. 

NID s 1 except when closed end and open 
end cases run separately. Then NID = 2. 

6 

KSKIP 

READS 

15 

= 0 
= 1 

For 1 case of COKPYN 

For successive cases using the same EOD 

output. 

11 

N4SCL 

READS 

15 


N4S0L = 0 when there are three solutions 
from EOD (one axi symmetric solution for 
each closed and open end cases and one 
crossflow solution for closed end). 

N4S0L = 1 when there are four solutions 
from EOD (solution for crossflow with 
open end is added). 

16 

NSPLT 

READS 

15 


Number of noise suppression devices 
(splitters): can equal zero. 

21 

KEYCOM 

READS 

15 

= 1 
* 2 
* 3 

Onbody data punch FLAG for data transfer 
Shroud Data 
Hub Data 

Shroud & Hub Data 

26 

KEYRK 

READS 

15 


The number of X's at which Vx, Vy, Vz 
is to be obtained. 

C-3 ONBODY INPUT CARD (4F10.5) 


COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

XRAG1 

READS 

F1C.5 


See Tape C-l. 

11 

XRAG2 

READS 

F10.5 


See Tape C-l. 
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TABLE C- A 


Use of Parameter XRAG 


1 

'EYCOM 

XRAG PARAMETERS 

MOTES 

1 

XRAG1 = Start value for shroud data 
(inside) 

XRAG2 = End value for shroud data 
(outside) 

XRAG3 and XRAG4 not used 

® If XRAG1 > 0 and 
and XRAG2 > 0 
— normal mode ana 
cnecKS for X valves on 
the shroud from inside 
to outside 



Q XRAGl < 0 
ABS (XRAGl) > XRAG2 
— Mode by which only X 
values on insiae surface 
of the shroud are 
obtained 



© XRAGl < i 
ABS (XRAGl) < XAAG2 
— Mode by wnich only X 
values on external 
surface of the shroud 
are obtained 


C-4 OFFBODY INPUT CARD (8F 10.5) (input only if KEYRK >0) 


1 

XOFFM(l) 

READS 

F10.5 

Offbody X at which Vx, Vy, Vz co be 
obtained. 

11 

XQFFM(2) 

READS 

F10.5 

Input KEYRK values 8 values per card. 

21 




71 




C-5 

NOISE SUPPRESSION SPLITTER 

INPUT (2015) (input only if NSPLT >0) 

COLUMN CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NSPB(l) 

READS 

15 

Number of right most point on 
Splitter. 

6 

NSPE(l) 

READS 

15 

Number of left most point on 
Spl itter. 
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11 

MSPB ( 2 ) 

READS 

15 

16 

tiSPE(2) 

READS 

15 


Input NSPLT values of NSPB and NSPE >8 values 
to a card. 

NOTE: The first Splitter is the one closest to the hub and tne last 

Splitter is the one closest to the shroua. 

C-6 FLOW CONTROL CARD#1 (8F10.2) 

ROUTINE 

COLUMN CODE FORMAT EXPLANATION 

Average axial velocity at the control 
station. Based on iive flow area, i.e., 
the flow area minus the area associated 
with the boundary layer displacement 
thickness. If VC = G.O the program will 
interpret this as a code and will 
calculate VC from WDOT. (To run a case 
with VC actually equal to zero set 
WDOT =0.0) (FT/SEC) 


11 

VINF 

READS 

F10.2 

Free Stream Velocity (FT/SEC) 

21 

ALFAF 

READS 

F10.2 

Angle of attack, 0.0 for free stream 
perpendicular to inlet axis (ap = X-90 u j 


degrees 

Total temperature, if PSTAT and TSTAT are 
read in (to be explained later), the program 
will calculate TTOTAL. If TTOTAL = u and 
PSTAT and TSTAT = 0, then TTOTAL = 516.67 
degrees k. 


41 

ELND 

READS 

F10.2 

Arbitary length used for normalizing 

(see KND card C-8 for further description; FT 

51 

YWING 

READS 

F10.2 

Not used in this program version YWING 
is the upper limit integration for surface 
forces. 

61 

UTIP 

READS 

F10.2 

Rotor tip speed (not required unless 
relative rotor inlet quantities are 
desired). 


31 TTOTAL READS 
F10.2 


1 VC READS 

F10.2 
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71 VA READS Sulk velocity at control station, i.e., 

F10.2 average inlet axidi velocity oased on 

geometric area. If VA = O.U, tne program 
will interpret this as a code and sec VA 
= VC. 

C-7 FLOW CONTROL CARD #2 (8F10.2) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

PT 

READS 

F10.2 

Total pressure, if PT = 0 program sets 
PT = 2116 Ibs/rt 2 ilos/FT 2 ) 

11 

PSTAT 

READS 

F10.2 

Static pressure, if input as zero, 
program will calculate (Ibs/FT 2 ). 

21 

TSTAT 

READS 

F10.2 

Static temperature, if input as zero, 
program will calculate (lbs/FT*). 

31 

WDOT 

READS 

F10.2 

Weight flow (input only if VC is 
input = 0.0). 

41 

DELQ 

READS 

F10.2 

Increment for value of stream function. 

C-8 GEOMETRY CONTROL CARD (2015) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NTHETA 

READS 

15 

Number of THETAS, where THETA is the 
circumferential coordinate. 

6 

NCLO 

READS 

15 

One rake must be chosen as tne control 
station. NCLO is the number of the first 
point on the rake. 

11 

NCHI 

READS 

15 

The number of the last point on the 
control station rake. 

16 

NX 

READS = 0 

15 ^0 

No XTEST values read. 

XTEST values must be input. 

21 

KND 

READS 

15 KND = -1, 

Flag used for non-dimensional i zing, 
ELND = YRISHR 


KND = 0, ELNO = 1.0 (no nondimensionalizingj 

KND = 1, ELND = YRISHR - YRIHUB 

KND = 2, ELND is the read in vaiue 

For the above the velocities are 
normalized by VC. 

KND = 4, same as KND = -1 but VA is used for 
normalizing velocities. 
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KND = 5 same as kwO = U, velocities normalized oy 
YA. 

KUO = 6 same as kND = 1, velocities normalized oy 
VA. 


C-9 THETA INPUT CARDS (8F10.2) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

THETA(l) 

READS 

F10.2 

Ci rcumferential coordinate in degrees. 
Note that the last theta read will De the 
data set used to calculate the Dounaary 
layer solution in VISCUS. 

11 

THETA (2) 




Input NTHETA values 

8 values per card 

C-10 X 

TEST CARDS 

(8F10.2) (NX 70) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

XTEST(l) 

READS 

F10.2 

Axial location of control surface, 
usually set equal to XRI 

11 

XTEST ( 2 ) 




Input NX 

values 8 values per card 

C-ll GEOMETRY DEFINITION CARD 

(8F10.2) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

XRi 

READS 

F10.2 

Value of X (axial coordinate) wnere 
surface distance is zero, usually equal 
to X at the control station. 

11 

YRIHUB 

READS 

F10.2 

i value on HUB at X = XkI 

21 

YRISHR 

READS 

F10.2 

Y value on SHROUD at X * XRI 

C-12 DATE CARD NO. 

, 1 (4X, 312, 

15) (Read only if KSKIP = G) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

5 

MO 

READS 

12 

Month of EOD run which produced 
potential flov solution. 
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7 

JA 

READS 

12 


Day of EOD run which produced 
potential flow solunon. 

9 

IR 

READS 

12 


Year of EOD run which produced 
potential flow solution. 

11 

NRUNND 

READS 

15 


Run number corresponding lo uate 
of computer program run aate given <ioove. 

C-13 

DATE CARD 

NO. 2 (4X, 312 

, 15) 

(Read only if KSKIP = 0 AND wib > ±) 

COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

5 

M02 

READS 

12 


See MO above 

7 

JA2 

READS 

12 


See JA above 

9 

IR2 

READS 

12 


See IR above 

11 

NRUNN02 

READS 

12 


See NRUNNO above 

* This card is input only if the open and closed potential flow solutions 
were run separately. 

The 

following cards are input for 

force and moment option. 

C-14 

INLET GEOMETRY CARD NO. 1 

(3A10) 

COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

TITLE 

READS 

(8A10) 


Problem identification card 

C-1S 

INLET GEOMETRY CARD NO. 2 

(12, 

flX, 12; 

COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

NINT 

READS 

12 


Total number of internal surface geometry 
points input 

11 

NOUT 

READS 

12 


Total number of external surface geometry 
points input 
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C-16 

INLET GEOMETRY CARD NO 

(2F10.0) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPlAnATIuk 

1 

XI(I) 

Rl ADS 
F10.0 

Internal surface geometry from ^CIkCc 
output (in combined routine) 

11 

YI(I) 

READS 

F10.0 

Data proceeds from XAAG* to nil itc 
location 

21 

OYDXIU) 

READS 

F10.0 



There will be NINT cards input (1=1, MINT) 

C-17 

INLET GEOMETRY CARD NO. 

4 (3F100) 

COLUMN 

COPE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

X0(I) 

READS 

F10.0 

External surface geometry from SCInCL 
output (in combined routine) 

11 

YOU) 

READS 

F10.0 

Data proceeds from hilite location to 
XRAG2 

31 

OYDXOU) 

READS 

F10.0 



There will be NOUT cards input (1=1, NOUT) 

NOTE: X = Axial Coordinate 

Y = Radial Coordinate 
OYDX = Surface Slope 

C-18 FLOW FIELD DA T A CARD NO. 1 (8A10) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

TITLE 

READS 

8A10 

Problem Identification 

C-19, 20 

FLOW 

FIEL DATA CARD NO. 

2 AND NO. 3 (NAMELIST) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

2-72 

SIN 

READS 

(NAMELIST) 

NAMELIST identified 


PAT4B 

TAMB 


Ambient Pressure, lbs/in t 
Ambient Temperature, °rt 
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VFLIGHT 

NACANS 

5-72 ACALPHA READS 

;.W£LIST) 


Free Stream Velocity, FT/jEC 
Nacelle Tilt Angle, oegraes 

Aircraft Angle of Attack, degrees 


SIDSLI? 


Sideslip Angle, degrees 


S* 

L* 

2 


Distance outboaru to moment reference 
point, inches 

Distance from high lite aiong inlet 
centerline to moment reference point, 
inches 

Terminator 


♦EXAMPLE: Reference point defined at the inlet 
high lite is S = 0.0, L = 0.0 

NOTE: 2IN must start in column 2. Flow conditions variaoles follow witn 

a comma separating one from another. If the data list is too large 
for one card, a second card can be used provided the data is listed 
after column 2. Do not end a card with a comma. Tne namelist is 
ended with a 2. 


C-21 FLOW CONDITION DATA CARO NO. 4 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

I 

INANG 

READS 

F10.0 

Angle between relative wind and inlet 
centerline 

21 

NPSIN 

READS 

12 

Number of circumferential positions von 
body) 

C-22 FLOW CONDITION DATA CARD NO. 

5 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NRAD 

READS 

12 

Number of circumferential positions 
ioff-body) 

11 

NCIRC 

READS 

12 

Number of data cards. This will equai 
the number of points per rake defined in 
SCIkCL 

21 

XFLOW 

REAOS 

P10.0 

Off -body axial station (race station) 
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C-23 FLOW COHO IT IOii DaTA CATO :.C . 0 


ROUTINE 

CCL'Jiif* CODE FORMAT 

1 RAD IUS ( I ) READS 

{aFiO.O} 


SPECIAL INSTRUCTIONS 


EXPLANATION 

Off body racial position. Tnese wili oe 
defined in the COhoYii output from one 
combined routine. Tnere will oe NCIkC 
number of RADIUS! I) data. 


The geometry data must cover the inlet boundaries defined in the combined 
input. These boundaries must be broad enough to include tne predicted 
stagnation point at each circumferential angle. If the boundaries are not 
broad enough, an error message will be printed. 

When defining the induced angle to calculate the forces and moments, tne 
aircraft centerline and the Nacelle centerline must be parallel. Therefore, 
in general, NACANG should equal 0.0 and ACALPHA should equal tne incidence 
angle. Then IHANG = ACALPAA. 
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i:4FUT FOR TkAWSG PRoGRrti 


This program takes data from tne COM3 YU program a.iu translates, scales, 
and sets up the aata for input to sub-program viscus. 

T-l TSAHSG CONTROL CFRD (F10.G) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 


1 

XBLT2 

TRANSG 

F10.0 

X station in tne iniet wnere tne 
layer solution ends. Usually it 
to XRI. 

t>o uiiGo ry 
is equal 
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liiPLT FOk VISCOUS ROUT I hit 


Since the program joes not use any constants wnicn depend on tne system 
of units (except the Free-stream total temperatures, any consistent set of 
units may be used. 

V-l SWITCH ARRAY CONTROL CARO (1515) 


COLUMN 

CODE 

ROUTINE 

FORHAT 


EXPLANATION 

1 

IST0(1) 

COTR 

15 

< 0 
§ C 

i 0 
= 999 

Data input on unit 51 
Data input on unit 5 
Input data listed 
Program returns at this point 

6 

ISTO ( 2 ) 

COTR 

15 

= -1 

Do not list profiles except at the 
X-location denoted oy ISTO ( 7 j or at the 
last X-location calculated if separation 
occurs. 




= 0 

Profiles are listed at each point to the 
profile end if U/V % 1.: if 0/V > i at 
end and WSTO(i) s 1 each succeeding 
point is listed unless IST0(4) p u then 
overshoot portion is suppressed entirely. 




= 1 

File lists only tne profile end values at 
the wall and boundary layer edge. 




= 2 

Profiles are listed on each point to the 
first point greater than W5T0 U / - 




= 3 

Profiles are listed on each point at tne 
UDIV interval to the profile end. 




> 3 

or LT < -1 only the value of the profile 
at the wall is listen. 

11 

ISTO (3) 

COTR 

15 


Used to compute the profile point list 
interval. List interval * ISTU(3j *JDIV 

16 

ISTC(4) 

COTR 

15 

* 0 

List for overshoot portion of profile 
(U/V) > 1 is suppressed. 

21 

IST0(5) 

COTR 

15 

= 1 

POP s 0 Lists momentum & energy equation 
boundary values for each iteration. 




= 0 

Oo not list transition statistics or 
viscosity functions at the wall. 

26 

ISTD(6) 

COTR 

15 

= 1 

Energy and momentum balances COFi , C0F2 
etc. together with their ratio are listed 
for all computed X-stations. 
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31 

ISTO ( 7 ) COT* 

> 

X-location numoer tnat a profile list is 


15 


required. 

36 

ISTO (8) COTR 

> 1 

Used to include a linear section of 


15 


X-station data from X = u to some X-iocation 


specifi ad by slope test. Tne end point of the 
linear section after the interpolation or new 
X points equals IST0(3< * X-station nuuoer of 
the last station input Oefoi-e tne 
interpolation. Tne total nuraaer of X-station 
points will be increased by ISTO {3; Inis 
option should oe used only when the si ope or 
the input data increases somewnere between a = 
C and the first maximum. 


41 

IS1D(S) 

COTR 

15 


Determines the end-point of the linear 
section. If the X-station data input > cnan 
ISTO ( 9 ? *2 begin the slope test. 

46 

ISTO (10) 

COTR 

15 


Not used by prig ram 


51 

ISTO(ll) 

COTR 

15 

* 1 
= 2 

One set of sta .ion data read in on tape. 
More than one ;et of station data read in 
tape 

on 

56 

ISTO (12) 

COTR 

15 


Non-dimension Oizes X-station data after 
computations if IST0(12) = 0 and JSTOiiQ) 
X-station data are non-oimentional ized oy 
last X-value. 

= 0 
the 

61 

IST0(13) 

COTR 

15 

= 1 

List full profile near zero shear 


65 

ISTO (14) 

COTR 

15 

= 1 

Turbulent restart occurs at point of 
laminar separation. 


71 

ISTO (15) 

COTR 

15 

= 0 

Station data input on tape 





t 0 

Station data input on cards. 


V-2 CONTROL CARD 

USTO(N0. 

1) (5E15.8) 


COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 


1 

HST0(1) 

COTR 

E15.8 


Limit for U/V usually = .SS99 


16 

WST0(2) 

COTR 

E15.8 


Used to limit amplitude of profiles 
U/V for calculation. 
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31 

WST0{3) 

COTR 

£15.8 

Used to recalculate koT(.> if 10? = *♦ 
tfSTuioj =1/ (i.) 

46 

WSTO (4 ) 

COTR 

£15.8 

Used only if [nX * > x 

61 

WST0(5) 

COTR 

E15.8 

Used as scale factor for X( I j jtfore 
computations. 

V-3 CONTROL CARO WSTOOJO. 

1) (5E15.8) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

U$TO(6) 

COTR 

E15.8 

See WSTD(4) 

16 

WSTO(7) 

COTR 

£15.8 

See WST0(4) 

31 

WST0{8) 

COTR 

£15.8 

Scale factor for U(I) 

46 

HST0(9) 

COTR 

£15.8 

Scale factor for GPiI) 

61 

WSTO(IO) 

COTR 

£15.8 

If NSTO(iG) = 0 and ISTOu^l = 0 set 
WSTO(IO) = X( IX ; i IX raters to last computed 
station. 




If WSTGClOj X 0 X is non-oiraensioned after 
computation and X(l) is scaled XU) = 

X( I )/WST0(iu) . 

V-4 CONTROL C/WD WSTO(NO. 

3) (5E15.8) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

WSTO(ll) 

COTR 

£15.8 

Used to compute RDTil) if IOP = 7 & DuP 
> 0 RDT(i) = U(I)0T(ii wnere NSTOim - v 
WSTOili) {if WSTO(li) = o 
program sets WSTO(ll) = .1564 x iu -3 } 

16 

WST0(12) 

COTR 

£15.8 

Not used 

31 

WST0(13) 

m o 
>-> o 

00 

Value of TURB(I) in the transition region 
for which turbulent restart will De allowed. 

46 

WST0(14) 

COTR 

£15.8 

= 0 Set PRT(turoulent prandil no) = 1.0 

40 PRT = WSTO ( 14 ) 
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61 WST0(15) COIR 

E15.8 


= 0 iMolecular pranotl no. = .78 
4 0 PR = WSTO l 15) 


V-5 CONTROL CARD W3T0(NG. 4} (5Ei5.3) 


COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

WSTO (16) 

COTR 

E15.8 


Conversion factor for U(I) 4 o. Then 
U(I) = WSTO (16) * Uili 

16 

WSTO (17) 

COTR 

E15.8 

4 0 

Than OBC scaled by WSTO i 17 1 

31 

WSTO (18) 

COTR 

E15.8 


Sonic velocity corresponding to static 
temperature at start of bounaary layer 
calculation. 

46 

WST0(19) 

COTR 

E15.8 


Static temperature at start of boundary 
layer calculation. 

61 

WSTO (20) 



Not used by this program. 

V-6 TITLE CARD (18A4) 



COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

LABEL 

COTR 

18A4 


Titie for this case 

V-7 C0N1R0L CARD NO. 1 FOR 

INPUT PROFILES (1515) 

COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

JOIV 

COTR 

15 


Number of suointervals of f'vn) and g'u 
input profiles, uDIV t 300 

V-8 CONTROL CARD 

NO. 2 FOR 

INPUT 

PROFILES (1515) 

COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

JY 

COTR 

15 


Total number of n values of f\n) and 
g'(n) input profiles, oY * 30u 

V-9 PROFILE INPUT CARDS (6F10.5) 


COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

YY(1) 

COTR 

F10.5 


First profile value. 
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11 

YY(IMUIV) 

COTR 

Fiu.5 

Second profile value. 

21 

YY( 1+2J0IV) 

COTR 
Flu. 5 

Third profile value 


etc. input until YY(JY) 

is reached. Input 6 values per card. 

V-10 

CONTROL CARD FOR F'iV) INPUT 11515) 

COLUMN 

l CODE 

ROUTINE 

FURMAT 

EXPLANATION 

1 

JEF 

COIR 

15 

Total number of f'(n) values or. the 
input profile. 

V-ll 

f'(V) INPUT CARDS (6F10.5) 


COLUMN CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

FP( 1) 

COTR 

F10.5 

1st f*(n) profile value 

11 

FP(1 + JUIV) 

COTR 

F10.5 

2nd f'(n) profile value 

21 

FP(1+2JDIV) 

COTR 

F10.5 

3rd f * (n) profile value 


etc. input until FP(JEF) 

is reached. Input 6 values per card. 

V-12 

CONTROL CARD FOR g' {*?) INPUT (1515) 

COLUMN CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

JEG 

COTR 

F10.5 

Total number of g'U) values on 
input profile. 

V-13 

q'(*>) INPUT CARDS (6F10.5) 


COLUMN CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

GP(1) 

COTR 
FIO. 5 

1st value of g'(n) profile value 

11 

GP(1 + JDIV) 

COTR 
FIO. 5 

2nd value of g'(n) profile value 
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21 GPU+2JDIV) 0073 
FIG. 5 


3rd value of g'i ) profile value 


etc. input until GPUEG) is reached. Input 6 values per card. 
V-14 PROGRAM OPTION C^D (1515) 

ROUTINE 

COLUMN CODE FORMAT EXPLANATION 


1 IGP COTR 
15 


6 MOP COTR 
15 


11 OOP COTR 
15 


Initialization option number: 

If IOP.EQ. 1,2,3 the input profile is the 
starting profile. 

If IOP.EQ.4,5,6,7, the starting profile is 
calculated using the input profile as a rougn 
guess. 

If I0P.EQ.1, the starting flow is laminar ana 
(dU/dx)^/^ is known. 

If I0P.EQ.2, the starting flow is turbulent 
and 6*(dp/dx) Ty is known. 

If IOP.EQ.3, the starting flow is turbulent 
and idU/dx)6*/U is known. 

If I0P.EQ.4, the starting flow is similar 
laminar. Omy on this option do calculation* 
start from the beginning of the ooundary layer 
growtn. 

If I0P.EQ.5, the starting flow is laminar ana 
(dU/dxU* 2 /*^ is known. 

If I0P.EQ.6, the starting flow is turbulent 
and 6*(dp/dx)/ry is known. 

If I0P.EQ.7, the starting flow is turbulent 
and (db/dx) 6*/U is known. 

Designated options on the g'(»?i profile. 

The option number is determined from the 
desired method of obtaining the g'vi?i profile 
whether the flow properties are variaole or 
not (see table III). 

Controls the interpretation of either the 
input velocity or the Mach number input aata. 
The option number is determined from caole IV. 
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ORIGINAL PAGE W 
OF POOR QUALITY 


Flow properties 

Method uf obtaining g ’(</ j protiie 


Variable 

-1 

1 

h° assumed constant a.ui equal to h® throughout the 
layer. g'(n) = 0.0. 

-2 

2 

Z' w = CBC = [h® - h w (x)]/(h® - h r ) is the wall 
boundary condition imposed on the energy equation. 

-3 

3 

g” = -(GDC)d,^/T w * -S tr d,/T w ts the wall 
boundary condition imposed on the energy equation 
where S, r = R L S, r = (*£ ' h r )] for 

laminar similarity starting solutions and 
5tr = S tr = *w/[Pe U ( h e * h r)] fcr 311 others - 


TABLE III. ALTERNATIVE VALUES OF MOP 


Interpretation of input U jn (I) 

Interpretation of input, 
fj n (tl) profile 

Mach number input, 

a U in (I) - M(Q 

Velocity input, 

\ a 'V(» 

-1 

1 

A i' la (n) - (U - ui/U 



V = V, 



w in » 

-2 

2 

\ a ll) ~ (P e V -pu)/(P e U) 





^Thesc Interchanges take place at the beginning of the program 
and thereafter Lf(I) ar.d f'(l) have their conventional mean- 
ings, whereas VU(I) represents (P w /P e )V w throughout the 
calculation for DOP=±2. 

TA3LE IV. ALTERNATIVE VALUES OF OOP 
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ORIGINAL PAGE jS 
OF POOR QUALITY 


16 

10 

COTR 

15 

Determines type of flow: 

IQ.EQ.-l for axisymmetric flow on insiae surface 
formed from radii drawn from buoy axis. 




I0.EQ.1 for axisymmetric flow on outsiae surface 
formed from radii drawn from body axis. 




IO.EQ.O or blank for planar flow. 

21 

TOP 

COTR 

15 

If TOP.cQ.l, transition is co oe 
calculated by the program. 

26 

POP 

COTR 

15 

Print option for listing data. 

If TOP.EQ.i, skip detailed list of output. 

V-15 

CASE INITIALIZATION CARD 

(8F10.5) 

COLUMN 

CODE 

RCUTINE 

FORMAT 

explanation 

1 

MR 

COTR 

F10.5 

Reference free-stream Macn numoer generally 
equal to the potential flow Mach number at start 
of boundary layer calculation: 




If IOP.EQ.4, HR=Ml2), and if iOP.NE.a, 
MR=M(1). 




If IOP.EQ.7 and MUi.EQ.O.u, i-ik is reset co 
0.001. 

11 

DT(1) 

COTR 

F10.5 

Displacement thickness at start of boundary 
layer calculation: 




If I0P.EQ.7 and DTlij.EQ.0.0, DTU)=0.0C± 
in the program. 

21 

RDT(l) 

COTR 

F10.5 

Reynolds number based on displacement thickness 
at the start of the boundary layer calculation: 



* 

If I0P.EQ.4, RDT(i) oecomes 

(xg-x^ )U(x 2 )/ and can be caiculiteu 

in the program (see WST0(3j input). 

31 

BS 

COTR 

F10.5 

Input for initial pressure gradient: 


If IOF.EQ.i, BS is (du/dxjo* 4 /!^. 
If I0P.EQ.2, BS is 6*idp/dxi/T^. 
If I0P.EQ.3, BS is idu/dx)6*/U. 
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6 




ORIGINAL PAG* » 

0F POOR QUALITY 


If ICP.EQ.4, 3S=i.O for rdUner-jKan 
stagnation point flow and oS=0.u tor 
Blasius flat plate f i ow . 

If IGP.Eg.B, BS is tub/dx; z* L /V rja . 

If IOP.Eg.6, BS is o*idp/ux j/T^. 

If 10P.EQ.7, BS is i du/ dx ) a */ u. 


41 

TO 

COTR 

F10.5 

Free-stream total temperature in 
degrees Kelvin. 

51 

BH 

COTR 

Enthalpy ratio, (h° - h r ;/n?. 

e e 

61 

FT 

COTR 

F10.5 

Free-stream turbulence, fraction of time 
during which the flow at a given position 
remains turbulent. 

Y-16 

STATION DATA CARD (8F10.5) 

{Read only if IST0I15) t 0; 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

X 

COTR 

F10.5 

Wall station locations (see Special 
Instructions for Preparing Input section;. 

11 

U 

COTR 

F10.5 

Free-stream velocity corresponding to 
each x-location. The arrays of X(I), ini), 
etc., are specified along the boundary surface 
for the downstream calculations. Tne .uaxiinum 
number of x-stations that can be computed is 
99. The last station data card snout u oe 
followed by an x-card with an x less than the 
previous x to switch out of the card read loop 

21 

TURB 

COTR 

F10.5 

Indicates fraction of flow that is 
turbulent. 


For laminar flow, TURB(I)=u.G; for turbulent 
flows, TURB ( I ) =1. 0. By changing TURBlI) from 
0.0 to 1.0, either abruptly or gradually over a 
distance of several x-stations, the effect of 
transition can be simulated. Also see W$7uii3) 


2-45 



ORIGINAL PAGE *9 
OF POOR QUALITY 


If the transition option is usea ii.e., 
TGP.=.i), a TllR-iii) input; is not needeu. 

If flow is all laminar, TURBu) input is 
not needed. 

If flow is turbulent, set TURo(I)=i.O at 
desired X(I; location. 

All other TURBII) at locations GT.XiI) vn 
be set to i.O by the program. 


31 

GBC 

COTR 
F10. 5 

Wall boundary condition on the energy 
equation, either g'(wall) or g"(wall). 

41 

RW 

COTR 

F10.5 

Radius of body surface in the same units 
as x. 




If 10=0, RWsO.O. 

51 

vw 

COTR 

F10.5 

Transpiration velocity in the same units 
as U. 

61 

sw 

COTR 

F10.5 

Uikuradse sandgrain rougnness scale in 
the same units as x. 

71 

cw 

COTR 

F10.5 

Longitudinal wall radius of curvature in 
the same units as x tif convex cw < 0; if 
concave cw > 0). 
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Special Instructions fcr Preparing Input 


Specification of profiles. - The input Intervals of n are subdivided by 
the input ofilY. For a laminar calculation the n step need net vary 
appreciably across the layer, and the product of JY and JDIV equal to ISO is 
usually adequate to define a profile. However, in a turbulent portion of the 
x-profiles a smaller step size should be prescribed close to the wall than is 
specified further out to Improve the accuracy in the “1 aw-of-the-wall " 
region. The n step size although variable from wall to boundary layer edge, 
remains fixed throughout the boundary layer calculations. The input step 
sizes should be specified for accurate results in both laminar and turbulent 
portions of the flow. The outer edge of the boundary layer in the n 
coordinate will not change appreciably as the calculations proceed downstream 
since n is normalized with «*, Smaller n spacing will be required throughout 
the layer if a very small x step is used. 

The specifications of the f'(n) and g'(n) profiles depends on whether 
these profiles are used as the starting profiles or the input profiles are 
recalculated. If the input profiles are to be used as specified, they should 
be compatible with the initial pressure gradient. For turbulent flow the 
profiles must be well defined in the "1 aw-of-the-wall" region. This region 
can often be specified by using same empirical "law of the wall*. If the 
input profiles are recalculated, the Initial profiles need not be very 
accurate since ihe calculated profiles usually converge very rapidly to their 
final value for almost any reasonable input profile. 

Specification of x-step size. - The sequence of the x- values defines the 
x-spacing at which calculations will be performed. The x-step size depends on 
tiie input station Mach number or velocity distribution. The x-step siz* 
should be inversely proportional to the magnitude of the slope of the velocity 
distribution. For large velocity gradients the. x-steps must be very small. 

The x-step size is most sensitive at the start of the boundary layer 

? rowth. However, if the slope of the U(X) as a function of x curve is linear, 
arger x-step sizes can be taken. The input options IST0(8) and !ST0(9) can 
be used to modify input x-station data to include a linear section at the 
start of the boundary layer growth. 

If the program is run separately, or if VAPE Is started in VISCOUS by the 
user, the x-step size and corresponding velocity profile can be specified by 
the user with the card iiput. Otherwise, subroutine TRAHS will automatically 
input the step size and velocity gradient as printed out at the end of COMBYN. 
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2.4 THREE DIMENSIONAL JET PROGRAM** 


HESS POTENTIAL FLO’J PRCGRAI! 


H-l TITLE CARD (PA1Q.G) 
ROUTINE 


COLIMN 

CODE 

FORMAT 


EXPLANATION 

I 

TITLE 

HINIT 

8A10 


Title of input case 

H-2 CASE CONTROL CARO (A5, 

8IS) 


COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

I 

CASE 

KINIT 

A5 


Case identifier 

6 

INOPT 

HINIT 


If INOPT >_ 0 set :ase = reads 



15 


Read jet data from previous run 
subroutine JETOLD 




= -I 

wooler jet method used 




--Z 

Heston jet method used 




<-Z 

Thames rectangular jet 

11 

I EXEC 

HINIT 

15 

=c 

=1 

Partial execution (used to verify-input) 
Full execution 




=-l 

Jet calculations only (not checked out) 

16 

MPR 

HINIT 

15 

=0 

=1 

No matrix print* 
Matrix printed 

21 

NGPS 

HINIT 

15 


Number of ground plane sections. 
Ground plane sections must be first 
non-lifting sections Input 

26 

JET 

HINIT 

15 


Number of independent jet systems 
(A jet system contains 1 or more jets) 

31 

L20 

HINIT 

15 

=0 

Panel data written on NT 20 except 
for sections with X7RPRN non blank 


>0 Panel data written on output unit 

<0 No panel data written on output file* 

regardless of XTRPRN flag 

** The control parameter normally used are denoted by an asterik (*). 
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36 

L21 

HINIT 

=0 

JET (wooler) print is or NT 21 



IS 





»1 

JET (wooler) print for final * 

iteration only is on output file 




>1 

All JET (wooler) print is on output file 

41 

L22 

KINIT 

=C 

Pressures and panel loads are written 



15 


on NT 22 except for sections with XTRPRN 
non-blank 




>0 

Pressures and panel loads written on 
output unit * 

46 

BL 

HINIT 

so 

No viscous solution 



15 

=1 

Viscous solution generated 

H-3 

CONTROL CARD FOR SECTION 

INPUT (1615) 



ROUTINE 



COLUMN CODE 

FORMAT 


EXPLANATION 

1 

NLSEC 

HINIT 


Number of non-lifting sections 



15 



6 

LIFSEC 

HINIT 


Number of lifting section S 12 



15 


11 

LASWAK 

HINIT 

>0 

Semi-infinite wake not used ) See 



15 


! Volume I 




=1 

Semi-infinite wake used* ) Section 4.C 

16 

NOFF 

HINIT 


Number of off body points 



15 



21 

IG 

HINIT 


Ignored panel flag (rcr panels which are 



15 


not exposed to the flow field) i.e. 
wing pane's under pylons, etc. 




=0 

*1 1 panel s used* 




-i 

Some panels ignored (Panels defined will 
be removed from solution matrix.) 

26 

MOf€NT 

HINIT 


Moment origin Input flag 



15 

*0 

No moment center input 




*1 

Moment center coordinates input 

31 

NSYMI 

HINIT 


X-Z symmetry flag 



15 

sO 

No symmetry 




=1 

Symmetry in X-Z planes 
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41 

INLETS 

HINIT 

15 

Number of inlets to be input 

4f 

I NIVEL 

HINIT 

15 

Inlet velocity input flag 
SO One inlet velocity per 
inlet input. This velocity 
used at all inlet control point. 




>0 Inlet velocities determined in 
inlet analysis routine at all inlet 
control points 



NOTE: NLSEC ♦ LIFSEC ♦ INLETS S 24 

H-4 REFERENCE DATA CARO (8F10.0) 


COUIttl 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

REFA 

HINIT 

F10.0 

Reference area for force coeff . 
(Normally wing area) 

11 

SPAN 

HINIT 

F10.0 

Reference span for moment coeff. 
(Normally wing span) 

21 

CBAR 

HINIT 

F10.0 

Reference length for moment coeff. 
(Normally wing mal) 

H-5 MACH NUMBER CARD (8F10.0) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

AM 

HINIT 

F10.0 

Input free stream mach no. 

H-6 ANGLE OF ATTACK CONTROL CARO (1615) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

IATACK 

HINIT 

15 

Number of angle of attack cases to be 
Input s 10 

H-7 ANGLE OF ATTACK CARDS (6F10.0) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

ALPHA(l) 

HINIT 

F10.0 

Angle of attack No. 1 
(Input In degrees) 

11 

PSI (1) 

HINIT 

F10.C 

Angl e of yaw No. 1 
(nose left is positive) 
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21 

ALPKA(2) 

KINIT 

F10.C 

Angle of attack No. 2 

31 

PSI (1) 

HINIT 

Angle of yaw No. 2 


etc 

F10.0 


Input IATACK valves of ALPHA and PSI 
H-8 KCffNT REFERENCE CARD (8F10.0) (INPUT ONLY IF HOMENT *0) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

ORIGNX 

HINIT 

F10.C 

X coordinate of moment center 

11 

ORIGNY 

HINIT 

F10.0 

Y coordinate of moment center 

21 

CRIGKZ 

HINIT 

F10.0 

Z coordinate of moment center 

H-9 CONTROL CARO FOR NON-LIFTING PANELS 

(1615) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NLS1R(1) 

HINIT 

17 

Number of strips in non-lifting 
section No. 1 

6 

NLSOR(l) 

HINIT 

15 

Number of panels per strip in 
nonlifting section Ko. l 

11 

NLS1R(2) 



16 

NLSCR(2) 




etc 




Input NLSEC values of NLSTR and NLSCR (NLSOR s 99) 

H-10 CONTROL CARO FOR LIFTING PANELS (1615) 

COLUMN 

COPE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

IXFLAG(l) 

HINIT 

Extra strip flag for lifting section 


15 No. 1 - This defines a strip of elements 

which are completely enclosed within 
another boty. No boundary conditions or 
sources are solved for on the elements of 
this strip. 
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s C no extra strip input 
* 1 first strip is extra 
a 2 first and last strips are extra 
a 3 last strip is extra 


6 

NSTRIP(l) 

HIHIT 

15 

Number of strips in lifting section 
No. 1 (includes extra strips) 

11 

NSORCE(l) 

HINIT 

15 

Number of source panels per strip in 
section No. 1 (does not include wake 
panels) 

16 

NUAKE (1) 

HINIT 

15 

Number of wake panels per strip in 
lifting section No. 1. 

21 

IXFLAG(2) 

HINIT 

15 


26 

NSTR1P(2) 

N 


31 

KSORCE (2) 

■ 


36 

HUAKE(2) 

m 



etc 

Input lifsec values of I3FLAG, MSTRIP, KSORCE and NHAKE 
NOTE: MSTRIP (I) S 20 and KSORCE (!) ♦ NUAKE (I) s 99 

H-ll CONTROL CARD FOR IGNORED PANELS (1615) (INPUT OMIY IF 16 >0 AMD LIFSEC >0) 
ROUTINE 

CCLWN CODE FORMAT EXPLANATION 

1 IG1(I,J) HINIT First ignored panel on lifting strip I 

15 

6 IGN(I,J) HINIT Last ignored panel on lifting strip I 

IS 

11 IGI(I.J) 

16 IGN(I.J) 
etc 

Repeat for Jal, LIFSEC, M, NSTRIP(J) 

H-12 CONTROL CARD FOR INLET 6E0«TRY (1615) 

ROUTINE 

COLUMN CODE FORWT EXPLANATION 

1 NLSTRK(K) HINIT Number of strips in inlet section K 

15 
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Hurter of panels per strip In ir.let 
section K 


£ ?;L£CR(K) HIMIT 

15 

Repeat for K * 1, Inlets 

H-13 GECIE7T.Y CCORCINATE IKFUT CARD NO. 1 (6F10.0, 2A10) 

This card reads in first card of each section with appropriate 
identifiers 


COLUMN 

COCE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

XB(1) 

INPUT* 

FIO.C 

X Coordinate of point 1 of given section 

11 

YB(1) 

INPUT 

F1C.0 

Y Coordinate of point 1 

21 

ZB(1) 

INPUT 

F10.0 

Z Coordinate of point 1 

31 

XB(2) 

INPUT 

F10.0 

X Coordinate of point 2 

41 

YBC2) 

INPUT 

F10.0 

Y Coordinate of point 2 

51 

ZB(2) 

INPUT 

F10.0 

Z Coordinate of point 2 

61 

SEC ID 

INPUT 

AlO 

Section identifier 

71 

XTRPRN 

INPUT 

AlO 

Print flag 

H-14 GEONETRY COORDINATE INPUT CARDS NO. 2 (6F10.0) 


These cards read in the 

rest of coordinates for each section 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

XB(I) 

INPUT* 

F10.0 

X Coordinate of point I 

11 

YB(I ) 

INPUT 

F10.0 

Y Coordinate of point I 

21 

ZB(I ) 

INPUT 

F10.0 

Z Coordinate of point I 
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NCTES: The gecnetry coordinate input card no. I is input at the first of each 
section followed by the geometry coordinate input cards no. 2. The 
nunber of no. 2 cards ir^ut depend on the number of panels and strips 
in each section. 

The gecnetry coordinate cards are input for each section for 
nonlifting, lifting, and inlet sections in the following required 
order: 

o All nonlifting sections 
o All lifting sections 
o All inlet sections 

♦The above gecnetry cards are read In subroutine LSTINP which Is 
called by subrouting HINIT. The data is read from unit 5 (input on 
cards, tape) in an 8A10.0 Format and transferred to unit 14 for use 
in the rest of the program, (subroutine input). 

H-15 OFF BODY GEOFCTRV INPUT (6F10.0) (Read only if NOFF >0) 

ROUTINE 


COLUMN 

CODE 

FORMAT 

EXPLANATION 

1 

XC(I) 

INPUT* 

F10.0 

X Coordinate of off-body pt 

11 

van 

INPUT 

F10.0 

Y Coordinate of off-body pts 

21 

ZC(I) 

INPUT 

F10.0 

Z Coordinate of off-body pts 

31 

xc(in) 

INPUT 

F10.0 

X, Y, Z Coordinates 

41 

YC(in) 



51 

ZC(I-l) 




Repeat until NCFF coordinates are input 


♦See Note Directly Above 


H-16 

INLET VELOCITY CARDS (6F10.0) 


(Input only if 

inlets >0 and INLVEL 

s 0) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

VINLET 

FCOM 

F10.0 

Inlet Velocity Ratio 
VINLET/ V„ 


Repeat 

this card for each angle of attack and each inlet 
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H — 1 7 INLET VELOCITY CARPS (SF10.3) 




ROUTINE 

COLUMN 

CCCE 

FORMAT 

1 

YCENT 

IhLINT 

F10.3 

11 

ZCENT 

INLINT 

F10.3 


(Input only if INLETS >0 and IKLVEL >0) 
EXPLANATION 

Y Location of Center of inlet circle 
Z Location of Center of inlet circle 
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WESTCN JET METHOD 


W-l JET CONTROL 

CARO (1615) 

(Read only If Jet >0 and INCPT = -2) 

COLUMN CODE 

ROUTINE 

FORMAT 

EXPLANATION 

i w 

INPUTW 

Control flag for a of jets 


15 

= 1 1 jet 

6 

INPUTW 

NJ s 2 2 jets 

Control flag for wake calculation 


15 

s 0 wake calculation used (only good 


for a = 0) 

s 1 wake calculation Ignored 


W-2 WESTON CONTROL CARD (8F10.3) (Read only If Jet > 0 & INOPT = -2) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

VINF 

INPUTW 

F10.3 

Free stream velocity, V^,, FT/SEC 

11 

XMAX 

INPUTW 

F10.3 

Max value of X/D in jet centerline 
and vortex path. (End of jet should 
be aft of configuration) 

21 

DX6T 

INPUTW 

F10.3 

Growth parameter used in jet 
centerline and vortex path. 

31 

DEXTAX 

INPUTW 

F10.3 

Initial stepsize of jet aX/D. 

41 

ZGP 

INPUTW 

F10.3 

Ground plane Z coordinate 

51 

GPI 

INPUTW 

i 0.0 

F10.3 * 0.0 

Control flag for STO calculation 
Wall jet calculated 
No wall jet calculations 

W-3 WESTON JET 

IDENTIFICATION CARD 

(Jet 1) (8F10.3) (Read only if Jet > C i 


INOPT « . 

■21 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

R 

INPUTW 

F10.3 

Jet to free stream velocity ratio 
Vj/V* . 

11 

PHI 

INPUTW 

F10.3 

Jet exhaust angle relative to 
Y-Z plane (Jet 1) 

21 

PSI 

INPUTW 

F10.3 

Jet exhaust angle relative to 
X-Z plane (Jet 1) 
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31 

XJ 

INPUTW 

F10.3 

X location of let 1 

41 

YJ 

INPUTW 

F10.3 

Y location of Je. 1 

51 

ZJ 

INPUTW 

F10.3 

Z location of Jet 1 

61 

DIAJ(l) 

INPUTW 

F10.3 

Diameter of jet exit Use consi start 
units) 

W-4 

WESTON JET IDENTIFICATION CARD (Jet 
-nmrrr — 9 — cm 

2) (8F1C.3) (Read only if Jet >C, 


COLUMN 

^ 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

RR (2) 

INPUTW 

F10.3 

See R above 


11 

PHI{2) 

INPUTW 

F10.3 

21 

PSI(2) 

INPUTW 

F10.3 

31 

XJ(2) 

INPUTW 

no. 3 

41 

YJ(2) 

INPUTW 

F10.3 

51 

ZJ { 2 ) 

INPUTW 

F10.3 


See PHI above 
See PSI above 
See XJ above 
See YJ above 
See ZJ above 


>A1 1 for Jet exit no. 


2 


51 DIAJ(2) INPUTW 
F10.3 


See DIAJ(l) above 


Note: Cards W-l, W-2, W-3 and W-4 must be input for each angle of attack and 

for each jet system i.e. number of s£ts of W-l through W-4 = Jet x IATACK. 

(Jet is from card H-2 and IATACK is from card H-6). 
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WCCLER JET METHOD 


WO-l WOOLER CONTROL CARDS NO.l (4110, 4A10) Read only if J>0 

and Inopt * - 1 

and IMets * 1 (1st Jet) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

MULT 

JET3IN 

no 

Number of jets in configuration 
mult = 1 on 2 

11 

IGEOM 

JET3IN 

110 

Not used by this version of code 

21 

NPS 

JET3IN 

no 

Number of integration intervals per 
jet segment 3 £ NPS £ 10 

31 

NOIT 

JET3IN 

no 

Nun&er of iterations to be performed 
on mutual interference velocities 
between jets 




Iter = 2 normally (See Ref. 4-2, pp. 
18-20) 

(Leave blank mult * 1) 

41-80 

SUBTITLE 

JET3IN 

4A10 

Additional title for this 
particular case 

W0-2 

WOGLER CONTROL CARD #2 (4110) Read only If J>0, 




IN OPT * -1 and 
NJet > 1 (2nd Jet) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

MULT 

JET3IN 

110 


11 

IGEOM 

JET3IN 

110 i 

See above 

> 

i 

21 

NPS 

JET3IN * 

no | 


31 

NOIT 

JET3IN > 

110 
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WO- 3 

WCCLER JET DATA CONTROL (110, 7F10.5) Read only if MULT > C 

COLUMN 

COCE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

N 

JET3IN 

110 

Number of intervals to be used in the 
numerical integration of the jet 
center line, N £ 100 

11 

GS 

JET3IN 

110 

Interval size to be used in the 
integration of the jet center line 
as fraction of leading jet diameter 

21 

SIGH 

JET3IN 

F10.5 

Yaw Angle Degrees 

31 

THETA 

JET3IN 

F10.5 

Pitch Angle Degrees 

41 

FEE 

JET3IN 

F10.5 

Roll Angle Degrees 

51 

YGP 

JET3IN 

F10.5 

Height Above Ground, FT 

WO-4 

WOGLER JET NO. 1 

DEFINITION (8F10.0) Read only if MULT > 0 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

XJl(l) 

JET3IN 

F10.0 

X location of jet 1 exit 

11 

YJKl) 

JET3IN 

F10.0 

Y location of jet 1 exit 

21 

ZJl(l) 

JET3IN 

F10.0 

Z location of jet 1 exit 

31 

PHID(l) 

JET3IN 

F10.0 

Jet 1 exhaust angle i deg 

41 

PSID(l) 

JET3IN 

F10.0 

Jet 1 exhaust angle *deg 

51 

OJET(l) 

JET3IN 

F10.5 

Jet 1 exit diameter 

61 

VELJI(l) 

JET3IN 

F10.0 

Freestream to jet exhaust velocity 
Ratio, Jet no. 1 

71 

FAN(l) 

JET3IN 

F10.0 

This option not checked out. 
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WO-5 WOOLER JET NO. 1 MIXING VALVES (8FI0.0) 




ROUTINE 


COLUMN 

CODE 

FORMAT 

EXPLANATION 

1 

A1 

JET3IN 

Ratio of effective core diameter to 



F10.0 

jet exit diameter for annular or 
vaned nozzles. (For uniform jet set 
= 1.0. 

11 

B1 

JET3IN 

Jet mixing parameter for annular or 



F10.0 

vaned nozzles, (see Ref. 4-2 for 
details) 

UO-6 WOCLER JET NO. 2 DEFINITION (8F10.G) (Read only If mult >1) 



ROUTINE 


COLUMN 

CODE 

FORMAT 

EXPLANATION 

1 

XJ 2(1 ) 

JET3IN 


11 

YJ 2( 1 ) 

F10.0 

See Wooler Jet 1 data card above for 

21 

ZJ2(1) 


definitions 

31 

PH ID (2) 



41 

PSID{2) 



51 

DJET2(1) 



61 

*U2(1) 



71 

*N(2) 



WO-7 WOCLER JET #2 MIXING 

VALVES (8F10.0) 



ROUTINE 


COLUMN 

CODE 

FORMAT 

EXPLANATION 

1 

A2 

JET3IN 

See A1 and B1 above for definitions 

11 

B2 

F10.0 


WO-8 WOOLER COALESCED JET 

INPUT (8F10.0) 



ROUTINE 


CGLUMN 

CODE 

FORMAT 

EXPLANATION 


1 DIARAT JET3IN Empirical factor for coalesced jet. 

F10.0 If Q < 20° DIARAT = 1 If Q > 20° 



Note: Cards wo-1 through wo-8 must be Input for each angle of attack and for 
each Jet system. I.e. Number of sets of wo-1 through wo-8 = Jet x IATACK. 

(Jet Is from card H-2 and IATACK Is from card H-6) 
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RECTANGULAR JET METHOD 


R-l RECTANGULAR JET CCN7RC4. CARD (1615) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

ICONFG 

RECJET 

15 

Nozzle configuration specification 
s 1 blunt 
s 2 streamwise 

6 

NSIKK 

RECJET 

15 

Number of source/sink lines used in 
model , 1 < NSINK < 5 

11 

SNEGV 

RECJET 

15 

Number of segments into which each 
vortex curve is divided, _< 30 

16 

NSEGS 

RECJET 

15 

Number of segments into which each 
source/sink line is divided, <_ 15 

R-2 RECTANGULAR JET D£FINI T ICN CARD (6F10.3) 

COLUMN 

CODE 

ROUTINE 
FORMAT ' 

EXPLANATION 

1 

R 

RECJET 

F10.3 

Jet-to-freestream velocity ratio 

11 

OELTAJ 

RECJET 

F10.3 

Jet injection angle, degrees 
('vt used) 

2i 

S 

RECJET 

F10.3 

Length downstream of jet exit that 
vortex and source/sink line are to 
extend, normally 15 to **0 jet 
diameters. 

31 

XJ 

RECJET 

F10.3 

X Location of jet exit 

41 

YJ 

RECJET 

F10.3 

Y Location of jet exit 

51 

ZJ 

RECJET 

F10.3 

Z Location of jet exit 

61 

PHI 

RECJET 

F10.3 

Jet exhaust angle relative to 
positive Z axis 

71 

PSI 

RECJET 

F10.3 

Jet exhaust angle relative to 
X axis 


2-61 



R-2A RECTANGULAR JET NOZZLE CARO (gFlO.3) 


CCLUTN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

CE 

RECJE" 

Jet equivalent diameter 
4(Area) 

Oe s v 

11 

AR 

RECJET 

F10.3 

Jet nczzle aspect ratio 
(length/width) 

R-3 SOURCE/SINK STRENGTH INPUT CARO (8F10.3) (Input if nsink >0) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

SIGMA(l) 

RECTJET 
F10 .3 

Strength of source/si r.k for 
source/sink line no. 1 

11 

SI GMA ( 2 ) 

RECTJET 

F10.3 

Strength line no. 2 

21 

SIGMA(3) 

m 

Strength line no. 3 

31 

SIGKA(4} 

m 

Strength line no. 4 

41 

SIff!A(5) 

(Input NSINK valves) 

Strength line no. 5 


Note: Cards R-l through R-3 must be input for each angle of attack and for 
each Jet systenm. i.e. Number of sets of cards R-l through R-3 « JET x 
IATACK. (JET is from card H-2 and 1ATACK is from H-6). 
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VISCOUS SOLUTION 


HV-1 VISCOUS FLOW CONTRa VARIABLES 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

RI 

HINIT 
(BSE TUP) 
Free Format 

Reynolds number per foot 
Ri * u* 

V 

v * kinematic viscosity 


UI 

HINIT 
(BSE TUP) 

Freestream velocity (ft/sec) 



Free Format 

If input as zero, program will use 
Mach number input to Hess program 

HV-2 VISCOUS INPUT CONTROL CARO 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

IDEFLT 

HINIT 
(BSE TUP) 

= 0 Data input below 
= 1 Defaults variables 



Free Format 

TO: 

LG1S = 1 
LG17 * 0 
LG20 * 0 
LG32 * 1 
TI * 518.6 7°R 
SWP = 0 
RL = 1/12 
KPRINT . 0 
ITPTFL x 0 
INP2FL s 0 




See next card for explanation of 
these values. 

HV-3 VISCOUS CONTRa VARIABLES (1) 

(Input only if IDEFLT = 0) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

LGI5 

HINIT 
(BSE TUP) 
Free Format 

Boundary layer method selection 
* 0 2 dimensional boundary layer 

equations used. 


= 1 3 dimensional Infinite yawed 

wing equation used. 
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LG17 

HINIT 
(BSETUP) 
Free Format 


LG20 

HINIT 
(BSETUP) 
Free Format 


LG32 

HINIT 
(BSETUP) 
Free Format 

HV-4 

VISCOUS CONTROL VARIABLES (2) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

1 

TI 

HINIT 
(BSETUP) 
Free Format 


SUP 

KINIT 
(BSETUP) 
Free Format 


RL 

HINIT 
(BSETUP) 
Free Format 

HV-5 

VISCOUS CONTROL VARIABLES (3) 

COLUMN 

COPE 

ROUTINE 

FORMAT 

1 

KPRINT 

HINIT 
(BSETUP) 
Free Format 


ITPTFL 

HINIT 


(BSE TUP) 
Free Format 


Transition intern" ttency flag 
s C Transition from laminar to 

turbulent flow is 
instantaneous. 

s 1 Transition from laminar to 

turbulent flow is gradual. 

Compressibility flag 
= C Incompressible flow equations 

sol ved. 

= 1 Compressible flow equations 

sol ved. 

Profile print flag 
s 0 All boundary layer velocity 

profiles will be presented. 

si Print output summary data 

only. (No velocity profiles 
will be presented.) 

(Input only if IDEFLT = 0) 


EXPLANATION 

Free stream temperature. 
Too 


Sweep angle in degrees 
(input only if L615 = 1) 


Geometry scale factor. 

The X-Y coordinates are multiplied 
by this parameter before calculations 
are started. 


(Input only if IDEFLT = 0) 

EXPLANATION 

s 0 Output data from BSETUP defin- 
ing boundary layer cases sup- 
pressed. 

si All data output from BSETUP. 

s 0 Output data from boundary layer 
routine suppressed. 

s 1 Output data from boundary 
layer routine printed. 
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IJJP2FL 

HINIT 
(BSE TCP) 
Free Format 

* 0 
= 1 

Input data to boundary layer 
routine not printed. 

Input data to boundary layer 
routine is printed. 

HV-6 

VISCOUS TRANSITION LOCATION CARO 



COL WIN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

LG16 

HINIT 
(BSE TUP) 
Free Format 

Boundary layer transition location flag. 
* 0 Transition point is input (sta- 

tion number). 




= 1 

Transition point will be 
calculated by the Michel method. 

HV-7 

VISCOUS 

TRANSITION LOCATION 

(Input only if LG16 = C) 

COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

NXT 

HINIT 

Station number at which flow becomes 


(oSETUP) turbulent. For all laminar flow, set 

Free Transition NXT greater than expected number of 

points. For all turbulent flow, set 
NXT * 2. Program will not work if NXT 
set to less than 2. 


Note: Input HV-6 and HV-7 (if necessary) for all lifting sections at each 

angle of attack if ISECT = number of lifting sections. Then input 
(ISECT) (IATACK) values. 


2-65 



3.0 VAPE INPUT FLOW CHARTS. 

This section contains a series of flow charts to assist the user in 
deciding what cards need to be input for the case in question. The various 
input cards are explained in Section 2 and card set numbers are assigned for 
each of the options applicable i.e. E-l refers to the primarily control card in 
the executive system section. Each option is presented separately with the 
first subsection (3.1) indicating which options will be exercised. No attempt 
is made in this section to show program execution logic or internal data 
transfer logic. 

3.1 MAIN VAPE PROGRAM. 

In the VAPE main program control routine (program vape), the input 
determines which options are to be executed. Since the inlet analysis program 
consists of a larger number of options, these are also selected in this main 
routine. 

Figure 4-1 presents the main program flow as determined by the input 
parameters I PROG and IRPG2. As can be seen in the figure IPROG determines the 
starting option and IPRG2 and to a certain extent IPASS determines the final 
option. For example, IPROG a 1 and IPRG2 * 7 indicates that the program will 
perform all options available through the Hess three dimensional program, but 
does not indicate anything about jet solutions, which will be referenced later. 
Therefore, if only a three-dimensional Hess solution is desired then IPROG = 7 
and IPRG2 = 7. Since the input data required is input within each option, this 
chart also indicates the input required and the order of the input. For 
example, if IPROG = 2 and IPRG2 = 7 then the data for the options, ECO, COMBYN, 
TRANSG, VISCUS, GE0H0D and Hess are required in that respective order. 

It should be noted that options 1 through 6 (SCIRCL through GE0M00) are all 
inlet analysis options* 

3.2 STOCKMAN INLET ANALYSIS SUB-PRCGPAM INPUT. 

The charts presented in this section are for options 1 throught 6 of the 
previous section (Figure 3-1). They are input in the order specified on the 
primary control card. The cards are identified by a card set number as given in 
section 2.0. 

The logic of these charts are indicated in Figure 3-2 by observing that 
card IA-6 in input in the data set only if ENREED is not equal to -1, 10, or 1. 
Figure 4-2 shows the logic of input to option number 1 (subroutine SCIRCL). 

Figure 3-3 shows the primary logic of the axisynmetric potential flow 
program data input sequence. The subroutine call to part 1 is where most of the 
data In input. Subsequently part 1 calls BASIC 1 and BASIC 2 for further 
input. Two very important items must be noted in Figure 3-3: 

(1) Card AP-1 has an input quantity NIN, which Is the flag to tell the 

program where the Input data Is orginating. If NIN a 5 then the data 
for the potential flow routine must be Input on cards. If NIN a 20 
then the data will be transferred from subroutine scircl and no data 
will be input on cards. The normal case is NIN a 20. 
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(2) If NIN = 5 then two sets of data must be input to the potential flow 
routine. One for an open ended body and one for a closed ended body. 
See Reference 3-1 for further information. 

The input logic to the combination option of the inlet analysis method is 
presented in Figure 3-4. The chart is sufficient to indicate the flow if input 
data to this program. 

The next option is the TRANS program which takes data output from the 
combination routine and translates scales and sets up in correct format for 
input to the viscous routine. 

Since most of the data is transferred internally, only one value on one 
card is required in this routine as input. Therefore, no logic diagram is 
presented for this option. 

The last option in the inlet analysis program is the viscous option where 
the boundary layer effects are computed. Again this data input is very straight 
forward and is explained adequate y by Figure 3-5. 

3.3 HESS THREE DUCNSIONAL POTENTIAL FLOW PROGRAM. 

All of the data required for Hess program including the jet data and 
viscous data is read fn subroutine HINIT and transferred to unit HT14 for use by 
the various subroutines. Since the data input is easier to understand when 
placed in the correct logic flow if the program the subroutine names where the 
data is required (and read from unit MT14) are used as reference. The input to 
this part of Vape requires some careful planning to make efficient use of the 
program. 

The majority of the expense in a Hess run is the matrix formation and 
solution. The solution is such that the incorporation of additional right hand 
sides to the problem, or additional angles of attach are second order as far as 
computation time is concerned. Therefore it is prudent to run all of the cases 
that are required at one time. There fore, if multiple angles of attack are 
input with repeat values, i.e. a =0, 5, 10, 0, 5, 10 etc. then the inlet 
velocity matrix of values and the jet solutions can be altered to give a 
complete range of cases. A typical run might consist of 9 angles of attack 0, 

5, 10, 0, 5, 10, 0, 5, 10 and three different jet velocity cases with associated 
inlet velocities. A sample input for such a case is presented in a later 
section. 

The logic flow chart shown in Figure 3-6 should be studied carefully to 
avoid any problems. 

Since the Input to the jet and viscous methods are relatively simple and 
are read from the Hess program, the charts for these routines are presented in 
this section as Figures 3-7 through 3-10. 
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IPROO «l| 


IPROO - S\ 


IPROO - 3| 


iproq - 4 


IPROO • 5 


IPROO - 6 


IPROO - 7 


MPROG - l 
CALL 
SC IF, CL 


MPROG - S 
CALL 
EOD 


MPROG - 3 
CALL 
COMBYN 


MPROG - 4 
CALL 
TRANSO 


MPROO - 5 
CALL 
VISCUC 


MPROO - 6 
CALL 
OEOMOD 


MPROO « 7 
CALL 
iiKsa 


X is \ no! 

IPRC2 - 1>H 


, IPRC2“l 
i.rPASs-i > 


is \ 

IPR02-3 


IS > S v NO 
IPRC2-4 >- J 


^ IS NO 

IPRC2-5 >— * 


* 3 \ HO 
IPRCP-6 >— I 


^ i: ^ 
IPASS-l 


n 

3 a 


Figure 3-1* Mnin Program Logic Flow 
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original page 19 

OF POOR QUALITY 


Report No: NADC-77119-30 



Figure 3-2. In^wt logic flow for Geometry Module of Stockman's 
Inlet Analysis Program 
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Figure 3-3* Inmt logic flow for Axisymmetric Potential Flow 
Module of Stockman's Inlet Analysis Program 
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Figure 3-U. Input logic flow for Combination Module of Stockman's 
Inlet Analysis Program 
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Note: If ISTO(l) .GE. 0 than cards V-2 

through V-15 are Input on cards 
and transferred to Unit 31 for 
future runs 


Figure 3-5. Input logic flow for Viscous Module of Stockman’s 
Inlet Analysis Program 
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FIGURE 3-6. INPUT LOGIC FLOW FOR HESS THREE DIMENSIONAL POTENTIAL FLOW PROGRAM 
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Figure 3-7- Input logic flow for Vaught Cooler 
Jet Module 



Figure 3..S. Input logic flow for 
Vought /Weston Jet 
Module 



Figure 3-9. input logic flow for 
Vought Rectangular 
Jet Module 
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4.0 SAMPLE f\SE INPUT 

4.2 SAMPLE CASE FOR INPUT TO INLET ANALYSIS ROUTINE 

This section contains tMo sample input cases. The first example is 
representative of the Input required to analyze an inlet configuration to 
determine the inlet separation boundaries. This case, therefore. Includes 
input for the geometry routine, the potential flow routine, the coabi nation 
routine, and the viscous routine. Note that since the data is transmitted 
internally from one rou : ne to the next, that only a minimum amount of Input 
is required. For example, only one card is Input for the potential flow 
routine. 

The first case is for the QCSEE Inlet using the following data: 

Ahi/Ath = 1*37 2:1 Ellipse Internal Lip Shape 

D Hl/°M = -905 a « 40 degrees 

V* SO NTS 


4 * *09— 

•3CSES frtL/ATrtsi. 37 


a 

?i3 - 

E3C 



3 


li 

1.7 J» 

X® 

- ? . 7C - 

19 

43.61» 


— mu a *.ts- 

2*1 ELLIPSE 3-M./C*=.90S 

e 

t-3 13.615 


4.9 

5.9- 

2.45 5.95 


X.O 


2.0 


-4.445- 




-4.»4S 


INPUT POP 


11m 61$ 
J-i 

2.4 

1.9 

*.<i 

?S. 

ll.alS 

2.4 

2.4 

36.9 

2.4 

GEOMETRY 

ROUTINE 

1-5 

Si 
- J.S 
1-001 

36.0 

13.0 

■m 

11.619 

11.619 

6.9 

6.7 1.792 

4.991 



2.0 

3.C 

4.991 

1.702 

4.9*1 



8.0 

$•66? 

9-9C.& 

ue 

a a 

2.C 

“i ? 

JJJ 




?.S02 

6.4*'S 

1.5 

4*9 

2.5*2 

6. of 5 
6. .55 
6i ,55 

11.61$ 
Jf • 9 

6m%5$ 



STOP 

20 












INPUT FOP POTENTIAL ROUTINE 

3 9 

OCSEC OUT 

til 699. 

4ML/6TK-1.57 

6 9 

.69 213 296 39 

69 


1.4 

2*66. 

-a— 

lS9. 

2966. 

-9»1 J 

961. 

966. . 

15.197' 

1.9 

•29 

a.t 9.9 

INPUT FOB 
COMBINATION 




»l. 


4M. 


RorrmiE 


& 5.619 

«5am — 4 

ENO 
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The second example Is representative of the input required to analyze an 
inlet configuration to determine the forces and moments acting on the inlet. 
This case, therefore, includes all of the input data contained in Case 1, with 
additional input data added for the force and moment program. 

This case is also fcr the QCSEE Inlet using the following data: 

A Hl/ A TH * 1*46 2:1 Ellipse Internal Lip Shape 

^Kl^M 8 *905 a s 70 degrees 


* .70 


V s 75 KTS 



9.. 127.62* -23. 5*1.3 l.i 3.0 

23*V1) *2$ 


*“ DIPOT FOP 
COKBDIATIOH 
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4.2 SAMPLE CASE FCR IV.PVJ TC HESS THREE DIKENSICKAL PROGRAM 

This is a sample input case for the configuration of Reference 4.1. 

This case is defined to read in single values of inlet velocity for each 
subcase to be considered. In this run, a subcase is defined to be one 
combination of angle of attack, one jet velocity ratio, and one inlet velocity. 
There are nine (9) subcases defined in this input. The flags are set to run a 
case with the data for a Vought/Weston jet node!. The input data for the jet 
model are not shown in this listing, but are shown later under the data 
presented in Section 4.3. The data for the jet models would appear directly 
before the 999. card shown in the following listing: 
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6.0 VAPE SUBROUTINES 


The VAPE program system contains over 190 subroutines. These subroutines 
can be divided Into seven primary areas as shown below: 

o VAPE Main System 

o Hess Lifting Program 

o Viscous Module 

o Inlet Analysis 

o Vought/Meston Jet Model 

o Wooler Jet Model 

o Thames Rectangular Jet Model 

The subroutines assigned to the VAPE main system are general routines, l.e. , 
read and write routines which are used throughout the VAPE system. In Section 
6.1, a complete listing of all subroutines, the main area the routine Is 
associated with, and a short explanation of the purpose of the subroutine is 
supplied. In Section 6.2, a series of flow charts are given showing the 
relationship of the various subroutines. 

6.1 SUBROUTINE DESCRIPTION 

All of the subroutines contained In the VAPE system are listed on the 
following pages In alphabetical order. The main system component that the 
routine is associated with is given along with a short explanation of the 
subroutine. 


SUBROUTINE 

MAIN SYSTEM COMPONENT 

DESCRIPTION 

ADAMS 

Wooler Jet Model 

Solves a system of N order 
differential equations by means of 
a fourth order ADAMS predictor/ 
corrector method. 

AFORM 

Hess Lifting Program 

Performs the DOT product of the 
source velocities and/or onset 
velocities at each control point 
with the normal vector at that 
point forms the A(I, J) matrix. 

AREAA 

Inlet Analysis 

Calculates the circular area 
between SHROUD and HUB. 

AVEV 

Inlet Analysis 

Primary calculation routine In 
subroutine COMBYN. Calculates 
areas and velocities. 

AXISA 

Inlet Analysis 

Compute axl symmetric velocity 
components VX, VY, VT, CP. 

BALANC 

Wooler Jet Method 

Establishes initial conditions for 


merged jet from momentum balance 
of 2 merging jets. 
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SUBROUTINE 

MAIN SYSTEM COMPONENT 

DESCRIPTION 

BALNCW 

Vought/Weston Jet Method 

Establishes Initial conditions for 
merged jet from momentum balance 
of 2 merging jets. 

BASIC I 

Inlet Analysis 

Reads In Input data and does 
Initial program calculations for 
uniform flow. 

BASIC 2 

Inlet Analysis 

Read In Input data and as 
Initial program calculations for 
non-uniform flow. 

BITEST 

Wooler Jet Model 

Tests for blockage of 2nd jet by 
1st jet and checks for Inter 
section of the two jets. 

BOUNDL 

Viscous Module 

Main program for Cebecl -Smith 
boundary 1<tyer code used In 
viscous solution. 

B SETUP 

Viscous Module 

Reads In boundary layer Input data 
and transforms Hess output data to 
boundary layer Input requirements 
for Cebecl -Sm th code. 

BUFSHR 

VAPE 

System routine to allow buffer - 
sharl ng . 

CALCRS 

Viscous Module 

Computes elements of a new right 
hand side for Hess from the 
boundary layer displacements and 
velocities. 

CALVEL 

Vought/Westor Jet Model 

Calls routines to calculate 
Induced velocities. 

CDELTA 

Vought/Wescon Jet Model 

Calculates trigonometric values of 
jet Injection angle. 

CFCAL 

Wooler Jet Model 

Computes direction cosines of 
local jet coordinate system. 

CFCAL 1 

Wooler Jet Model 

Computes directional cosines at 
local jet coordinate system. 

COLSOL 

Hess Lifting Program 

Matrix solution routine. 

COMBYN 

Inlet Analysis 

Main control routine to calculate 
velocities for Inlet analysis 
routine with viscous effects. 

comow 

Hess Lifting Program 

Routine to control velocity 
solution after matrix is solved. 
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SUBROUTINE 

WIN SYSTEM COMPONENT 

DESCRIPTION 

CCH> 

Wooler Jet Model 

Computes V/V effective and tests 
for intersection of centerlines. 

CONST 

Inlet Analysis 

Calculates constants used in 
subroutine CONBYN 

COTR 

Inlet Analysis 

Reads In all input data for 
viscous routine. 

CROSS 

Inlet Analysi s 

Computes cross flow velocity 
components. 

CRPROO 

Vought/Ueston and 
Rectangular Jet Method 

Calculates cross product of two 
vectors. 

CUBIC 

Inlet Analysis 

Fits a cubic between 2 straight 
lines. 

GORVFM 

Meston Jet Model 

Forms jet or vortex curves. 

DBCQDU 

Inlet Analysis 

Computes the double integral of a 
table using splines. 

DBCQDU 

Wooler Jet Method 

Integration of F(X) from A to B 
using cautious adaptive ronberg 
extrapolation. 

DCOSJ 

Weston Jet Model and 
Rectangular Jet Model 

Computes direction cosines of 
local coordinate system. 

DCCSJI 

Rectangular Jet Model 

Computes direction cosines of jet 
deflection (routine identical to 
DCOSJ). 

DCSQDU 

Inlet Analysis 

Integrates a cubic spline between 
A and B. 

OERIV 

Wooler Jet Model 

Computes derivatives for ADAMS 
predictor 1 corrector method. 

DIRCOS 

Vought/Weston Jet Model 

Computes direction cosines between 
2 points. 

EDVS 

Viscous Nodule 

Calculates eddy viscosity for 
CEBECI-SMITH boundary layer method 

EINF 

Viscous Nodule 

Calculates transfonaed Y-grid 
points for Cebeci -Smith boundary 
layer method. 

ELIP 

Inlet Analysis 

Hastings approximation for 
elliptic integrals. 


C-2~ 
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SUBROUTINE 


.'AIN SYSTEM COMPONENT 


DESCRIPTION 


ENGY 

Viscous ‘todule 

Calculates solution to energy 
equation for Cebeci -Smith boundary 
layer method. 

EDO 

Inlet Analysis 

Hess axl symmetric potential flow 
main contrci program. 

ERF 

Weston Jet Model 

Error function calculation routine. 

EXCRCS 

Inlet Analysis 

Compute extra cross flow 
velocities for axi symmetric 
potential flow. 

FILE 

Inlet Analysis 

Lists the momentum and energy 
profiles on Y and the wall station 
parameters at specified X stations. 

FINOYS 

Inlet Analysis 

Finds value of Y for hub and 
s! x>ud at rake location - 
(Determines min. and max. height 
at rake). 

FIX 

Hooler Jet Model 

Limits mutually induced velocities 
to a maximum valve. 

FLP2 

Viscous Module 

Determines fluid properties for 
Cebeci-Smith boundary layer code. 

FMCAR 

Inlet Analysis 

Control program for force and 
moment calculations. 

FNCHON 

Inlet Analysis 

Compute moment force on Internal 
surface. 

FMDATN 

Inlet Analysis 

Reads data and Initializes 
calculation. 

FtfLUX 

Inlet Analysis 

Computes momentum flux. 

FMINP 

Inlet Analysis 

Reads data for force and moment 
calculations. 

FMINT 

Inlet Analysis 

Compute inside P-A integral In X 
direction. 

FMOUT 

Inlet Analysis 

Output routine for force and 
moment calculations. 

FNSTRH 

Inlet Analysis 

Find straight segment on hub and 
shroud. 

FONISO 

Inlet Analysis 

Find N for the supercircle 
F(N) = (X/a) n ♦ (Y/B) n -1 
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SUBROUTINE 

WIN SYSTEM COMPONENT 

DESCRIPTION 

FCRKCM 

Inlet Analysis 

Main control progran for force and 
nonent calculation. 

FRSTSH 

Inlet Analysis 

Find first straight segnent on 
shroud. 

GEOHOD 

Inlet Analysis 

Modifies ii.let geometry to account 
for boundary layer displacenent 
thickness. 

GETABC 

Inlet Analysis 

Conputes VI, V2, V3 and A, B, C 
frou input paraneters. 

GETMT 

Inlet Analysis 

Opens files. 

GROUNDH 

Inlet Analysis 

Calculates ground height. 

NEADK 

Viscous Nodule 

Prints header on each page of 
Cebeci-Smith boundary layer output. 

HESS 

Hess Lifting Program 

Main control progran for Hess. 

HINIT 

Hess Lifting Progran 

Reads all input data for Hess, 
including jet data. Lists it out 
and transfers it to Unit 14. 

ICSEVU 

Inlet Analysis 

Evaluation of a cubic spline. 

ICSICU 

Inlet Analysis 

Interpolator approxlnation by 
cubic splines with arbitrary 
second derivative end conditions. 

ICSSU 

Inlet Analysis 

Cubic spline data snoothing. 

INBLOCK 

Heston Jet Model 

Determination of equivalent 
velocity due to jet blockage at 
2nd jet exit. 

INCRNT 

Rectangular Jet Model 

Calculates increnent length for 
jet curves. 

INDEX 

Inlet Analysis 

A functional subroutine which 
compares or tabular function X 
with a functional Z to find the 
index of X such that X is the 
least tabular value of X greater 
than Z. 

INIR 

Inlet Analysis 

Determines Initial conditions for 
boundary layer calculations. 

INITMT 

VAPE 

Initial routine in VAPE - se** 
program c ontnoiAumainet* 1 ^ 
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SUBROUTINE 

MAIN SYSTEM COMPONENT 

DESCRIPTION 

INLIKT 

Hess Lifting Program 

Finds inlet CP's by interpolation. 

INPUT 

Hess Lifting Program 

Reads data from Unit 14 for Hess 
3-D potential flow routine. 

INPUT* 

Vought/Vieston Jet Model 

Reads input dat? for Heston 
routine. Also « eteraines jet 
spacing in terms of jet diameters. 

we 

Viscous Module 

Processes all input to the 
Cebecl-Smith boundary layer code. 

INS2 

Viscous Module 

Linear or quadratic interpolation 
routine. 

DiTEG 

Inlet Analysis 

Integrates sum by trapezoidal rale. 

INTEGA 

Wooler Jet Model 

Integration of the equations of 
notion for the jet path. 

INTPOL 

Inlet Analysis 

Interpolation and differentiation. 

INTSEC 

Vought/Ueston Jet Model 

Determines if merging of jets 
occurs and what direction cosines 
of jets 1 and 2 are at mergi ng 
point. 

ITFR 

Inlet Analysis 

Generates new initial FP and GP 
profiles for viscous solution. 

ITSR 

Inlet Analysis 

Main subroutine for computation of 
momentum and energy profiles. 

ITS2 

Viscous Module 

Integration routine in 
Cebecl-Smith boundary layer code. 

IVPF 

Viscous Nodule 

Generates initial velocity profile 
for Cebeci-Smith boundary layer 
solution. 

JDRCOS 

Vought/Ueston Jet Model 

Calculates direction cosines of 
freestream and jets. 

JETCET 

Vought/Ueston Jet Model 

Calculates jet centerline location. 

JETCT2 

Vought/Ueston Jet Model 

Jet centerline location for wake 
effects. 

JETOLD 

wooler Jet Model 

This routine reads jet induced 
velocities computed in an earlier 
run. 

JET 3 

Uooler Jet Model 

Main program for Wooler jet code. 
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SUBROUTINE 

MAIN SYSTEM COMPONENT 

DESCRIPTION 

JET3IN 

Wooler Jet Model 

Input routine and parameter 
initilization for Wooler jet code. 

JET3IT 

Wooler Jet Model 

Does Initial calculations, 
including jet blockage and 
Intersectional. 

JET3ZZ 

Wooler Jet Model 

Calculates jet induced velocities. 

LEM 

Inlet Analysis 

Calculates points on a leoiscate. 

LIFT 

Hess Lifting Program 

Calculates geometric quantities 
for elements in the lifting 
section. 

LIN1P3 

Rectangular Jet Model 

Lineary interpolation of 3 
independent variables. 

LSEP 

Viscous Nodule 

Determines when a laminar 
separation point, as predicted, 
will be a transition location or 
an actual separation point. 

LSTINP 

VAPE 

Lists all input cards. 

MATRIX 

Inlet Analysis 

Cooputes matrix for potential flow 
solution. 

NATSOL 

Hess Lifting Program 

Controls matrix solution routine 

(1) Col sol -Basic matrix solution, 

(2) Solmor- Solve for additional 
right hand sides with same matrix. 

MERJET 

Vought/Weston Jet Model 

Determines jet centerline of 
merged jets. 

MIRROR 

Inlet Analysis 

Mirrors the hub to obtain points 
on shroud (2-D inlets) 

MISNA2 

Inlet Analysis 

Matrix solutlon-Sledel Integra 
tlon method. 

MISI 

Inlet Analysis 

Matrix solution-direct Inverse 
method. 

MODIFY 

Wooler Jet Model 

Computes mutually induced 
velocities between jets. 

MOMX 

Viscous Module 

Finds solution of X-momentum 
equation in Cebecl-Smlth boundary 
layer code. 

MOMZ 

Viscous Module 

Finds solution of Z-momentum 
equation in Cebecl-Smith boundary 
layer code. 
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SUBROUTINE 


WIN SYSTEM CCMPGNENT 


DESCRIPTION 


NUINT 

Hooler Jet Model 

Computes direction cosines of 
modified freestream and modified 
free stream to freestream velocity 
rati o. 

NOIR 

Inlet Analysis 

Print boundary layer results 
summari zed. 

NEAR 

Hess Lifting Program 

Calculates source and dipole 
velocities of a lifting element in 
the near field- influe nee 
velocities for elements. 

NOEPTS 

Inlet Analysis 

Determines end points of rakes 
near hubs shroud and/or splitters. 

NOLIFT 

Hess Lifting Program 

Computes geometric quantities of 
the now lifting elements. 

OFFBDY 

Inlet Analysis 

Calculates off-body variables. 

OLCMFLO 

Hess Lifting Program 

Combines flow solutions and 
determines final velocities and 
pressures. 

ONBODY 

Inlet Analysis 

Calculates on-body variables. 

OTPT 

Viscous Module 

Outputs the results of the 
Cebeci -Smith boundary layer code. 

OUTPT 

Hooler Jet Model 

Transforms local coordinates to 
program coordinates. 

OUTPTI 

Hooler Jet Model 

Transforms local coordinates to 
program coordinates for segmented 
jets. 

PART 1 

Inlet Analysis 

Control for basic data and matrix 
formulation. 

PART 2 

Inlet Analysis 

Computer source density sigma by 
siedel iteration. 

PART 4 

Inlet Analysis 

Computes velocity components and 
print. 

PKUTTA 

Hess Lifting Program 

Computes vortex strength by kurta 
condition. 

PLANE 

Hooler Jet Model 

Computes intersection of plane 
with a line. 

PRELPS 

Inlet Analysis 

Calculates geometric propertlties 
in geometry portion of inlet 
analysl s. 
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DESCRIPTION 


SUBROUTINE 

."AIN SYSTEM COMPONENT 

PREP 

Inlet Analysis 

PRINT 

Hess Lifting Program 

PROFYL 

Inlet Analysis 

PRTOUT 

Wcoler Jet Model 

READS 

Inlet Analysis 

READ 1 

VAPE 

READ 3 

VA PE 

READ 4 

Inlet Analysis 

RECJET 

Rectangular Jet Model 

REGEN 

Inlet Analysis 

RETMT 

VAPE in general 

ROTATE 

Wooler Jet Model 

SCIRCL 

Inlet Analysis 

SDIST 

Vought/Weston Jet Model 
Inlet Analysis 

SEGMNT 

Wooler Jet Model 

SETMT 

VAPE 

SHFT 

Vlscus Module 

SIMQ 

Inlet Analysis 

SINK 

Vought/Weston Jet Model 


Set up tapes for MATSCL solution. 

Prints final output. 

Computes F(ETA) and G(ETA) in 
viscous program. 

Prints ut jet centerline data and 
induced velocities at all control 
poi nts. 

Reads all input. 

Reads block data from a unit. 

Reads multi block data from a unit. 

Reads data for force and moment 
calculations. 

Main program for rectangular jet 
method. 

Transfers data from Unit MT20 to 
MT29; part of viscous inviscid 
i nteraction. 

Tape manipulation. 

Rotates fixed coordinates to 
rotated or rotated to fixed. 

Main geometry routine for inlet 
analysis. 

Calculates surface distances. 


Establishes initial conditions for 
new segment of jet, for continu- 
ation, or a jet and Integrates 
jets. 

Sets up tapes. 

Provides Initial guess for each 
station in Cebecl -Smith boundary 
layer code. 

Solves set of simultaneous linear 
equations. 

Calculates Induced velocities to 
due sinks along jet centerline. 
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SUBROUTINE 

MAIN SYSTEM COMPONENT 

DESCRIPTION 

sintp 

Inlet Analysis 

Calculates geometric parameters 
used in combi nation routine. 

SLN6 

Viscous Module 

Determines surface distance from 
X-Y input for Cebeci -Smith boundary 
layer code. 

SLOPE 

Viscous Moaule 

Determines derivative DY/DX from X 
vs Y input in Cebeci-Snith boundary 
layer code. 

SOLMOR 

Hess Lifting Program 

Matrix solution routine. Solves 
basic Hess matrix with additional 
right hand sides. 

SOLVIT 

Inlet Analysis 

Matrix solution routine. 

SPL INS 

Viscous Module 

Interpolation by cubic spline. 

SRTNE 

Inlet Analysis 

Calculates surface distance used 
in combination routine. 

STRAIT 

Inlet Analysis 

Calculates straight line segments 
in SCIRCL 

STRLINE 

Inlet Analysis 

Calculates straight line segment. 

STRM. 

Inlet Analysis 

This routine calculates streamlines 
used in combination routine. 

SUPERC 

Inlet Analysis 

Geometry routine in SCIRCL routine. 

SUPRD 

Inlet Analysis 

Determines geometry of noise 
stqipression splitter plates. 

TEST 

Inlet Analysis 

Geometry routine used in SCIRCL 
geometry routine. 

TRANS 

Inlet Analysis 

Transition location prediction 
routine. 

TRANSA 

Wooler Jet Model 

Coordinate transformation routine. 

TRANSF 

Vought/Weston Jet Model 

Coordinate transformation routine. 

TRANSG 

Inlet Analysis 

Sets up data for viscous solution. 

TRANS3 

Wooler Jet Model 

Coordinate transformation routine. 

TRICK 

Inlet Analysis 

Smoothing routine. 

TRNSF1 

Rectangular Jet Model 

Coordinate transformation routine. 
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SUBROUTINE 

MAIN SYSTEM COMPONENT 

DESCRIPTION 

TRN2 

Viscous Module 

Computes location of boundary layer 
transition for Cebeci -smith code. 

UERTST 

Inlet Analysis 

Message generation. 

VAPE 

VAPE 

Main routine for VAPE program. 
Controls basic logic of program. 

VAROFF 

Inlet Analysis 

Calculates fluid properties off the 
body. 

VBARIT 

Inlet Analysis 

Solves VBAR combination Iteratively. 
Used In combination routine. 

VCOM 

Hess Lifting Program 

Velocities of configuration 
determined. Jet and inlet 
velocities determined through 
control of this routine. 

VELQC 

Wooler Jet Model 

Determines jet Induced velocities. 

VEL 1 

Wooler Jet Model 

Computes effective velocity ratio 
for downstream jet at exit. 

VFM.FT 

Hess Lifting Program 

Computes the induced velocities in a 
lifting strip. 

VFM1LF 

Hess Lifting Program 

Computes the induced velocities In a 
non-lifting strip. 

VFORM 

Hess Lifting Program 

Controls calculation of velocity on 
a body* plus forms v-matrlx. 

VIS 

Inlet Analysis 

Computes effective viscosity and 
conductivity. 

VISCOS 

Inlet Analysis 

Main control program for viscous 
routines. 

VORCUR 

Vought/Weston Jet Model 

Definition of vortex curves. 

VORT 

Vought/Weston Jet Model 

Calculate Induced velocities due to 
vortlcies. 

VORTVEL 

Wooler Jet Model 

Vortex Induced velocities 
calculated. 

WALJET 

Vought/Weston Jet Method 
(STOL) 

Calculates effect of wall -jet on 
ground plane. 

WESTON 

Vought/Weston Jet Method 

Main control program for Weston 
jet method. 
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SUBROUTINE 

MAIN SYSTEM COMPONENT 

DESCRIPTION 

UN EAR 

Hess Lifting Program 

Calculates dipole velocities of 
wake element In near field. 

WPUNCH 

Inlet Analysis 

Calculates geometric data and 
places on Unit 1 (MT20). 

URITE 1 

VAPE 

Write block data. 

WRITE 2 

VAPE 

Urite block data. 

URITE 3 

VAPE 

Write block data. 

URITE 4 

VAPE 

Urite block data. 

WRTXY 

VAPE 

Write data on Unit 20. 

XPROO 

Wooler Jet Model 

Confutes cross product of two 
vectors. 

XYZ 

Inlet Analysis 

Control for X, Y, Z matrices 
computation. 

XYZ1 

Inlet Analysis 

Compute X, Y, Z matrices for JJ 
less than or equal .08. 

XYZZ 

Inlet Analysis 

Compute X, Y, Z matrices using 
Simpson rule integration. 

6.2 SUBROUTINE FLOW CHARTS 



Flowcharts of the VAPE code showing relationships between subroutines 
are presented In this section. Figure 6.1 shows the main flow of the VAPE 
system. Figure 6.2 shows the logic for rutlne GEOMOD. Figure 6.3 shows the 
flow for the axlsymmetric Hess program. The viscous logic flow Is contained 
In Figure 6.4. The Input geometry formulation for the in’ :t analysis is 
contained In Figure 6.5. The routine which combines the flow fields In the 
Inlet analysis routine Is presented In Figure 6.6. The Hess routines and 
various jet methods are presented In Figures 6.7 through 6.9. In these 
diagrams, the subroutines are called from the other routines as shown. For 
presentation purposes, some of the reads, writes, and tape manipulation 
routines have been shown In the same block 
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FIGURE 6-1. FLOW 
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FIGURE 6-2. FLOW CHART OF GEOMOD SUBROUTINES 
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FIGURE 6-3. FLOW CHART OF EOD (AXISYMMETRIC HESS) SUBROUTINES 
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FIGURE 6-4. FLOW CHART OF VISCUS SUBROUTINES 


REPORT NO: NADC-77119-30 

















FIGURE 6-5. FLOW CHART OF SCIRCL SUBROUTINES 
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REPORT NO: NADC-77119-30 
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REPORT NO.: NADC-77119-30 
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FIGURE 6-8. FLOW CHART OF WESTON SUBROUTINES 
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REPORT NO.: NADC-77119-30 
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FIGURE 6-10. FLOW CHART OF VISCOUS MODULE 



















APPENDIX A - DESCRIPTION OF GEOMETRY PROGRAM FOR DEVELOPING HESS THREE 
DIMENSIONAL INPUT 

IN-1 PRIMARY CONTROL CARD (2A1, AS, 7A10) 

This card controls the main flow of the program. 


ROUTINE 


comm 

CODE 

FORMAT 

EXPLANATION 

l 

INST(l) 

VINPEX 

A1 

Selects module to be performed 
3 F Wing module 

3 j Junction module 

3 C Conic or circle module 

depending on INST(2) 

3 & Hesplot module - plots 

3 S Stop 

2 

INST(Z) 

VINPEX 

A1 

3 0 Conic routine 

3 I Circle routine 

3-80 

INST(3)- 

INST(IO) 

VINPEX 
A8, 7A10 

used to complete description of 
module desired not used in program 
except for output (i.e. foil, Juncpan 
Hesplot, Stop, etc.). 
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INPUT TO WING GENERATION ROUTINE 


This 

routine creates a wing panelled for direct input to the HESS program. 

F-L MAIN CONTROL OtKD U6I5) 


- 

COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

NPTS 

Foil 

15 

Number of coordinates on each surface 
of the airfoil 

6 

NFOILI 

Foil 

15 

= 0 
* 0 

NACA airfoil calculated for 
Inboard profile shape 
Airfoil coordinates for inboard 
profile will be input. 

11 

NFOILO 

Foil 

15 

= 0 
# 0 

NACA airfoil calculated for 
outboard profile shape 
Airfoil coordinates for out- 
board profile will be input. 

16 

IND 

Foil 

15 

= 0 
* 0 

Divide wing into NS equal 
panels in spanwise direction 
Read in non-dimensional span- 
wise values (N = y/b/2) 

21 

NS 

Foil 

15 

Number 

of spanwise airfoils to be used. 

26 

IDIHED 

Foil 

15 

» 0 
* 0 

No dihedral angle considered 
Dihedral angle will be used in 
defining wing 

31 

ITWIST 

Foil 

15 

» 0 
* o 

No wing twist considered 
Wing will be twisted 

36 

IVERT 

Foil 

15 

- 0 
* 0 

Horizontal tail, wing, canard 
will be defined 
Vertical tail or pylon will be 
defined 

P-2 WING DEFINITION 

CARD (7F10.0) 



COLWH 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 

1 

CROOT 

Foil 

F10.0 

Root chord 

11 

CTIP 

Foil 

F10.0 

Tip chord 

21 

SWEEP 

Foil 

F10.0 

Leading edge sweep angle deg. 



COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 


31 

B0V2 

Foil 

F10.0 

Distance from root to tip (b/2) 


41 

(TRYT.E 

Foil 

F10.0 

X 

Location of root leading 

edge 

51 

CHYLE 

Foil 

F10.0 

Y 

Location of root leading 

edge 

61 

rp?Tj; 

Foil 

F10.0 

Z 

Location of root leading 

edge 

F-3 DIHEDRAL DEFINITION CARD (7F10.0) (Read 

only if IDIHEE>0) 


comm 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 


1 

DANGLE 

Foil 

F10.0 

Dihedral angle in degrees 


11 

YD 

Foil 

110.0 

Y 

location of rotation for dihedral 

21 

ZD 

Foil 

F10.0 

Z 

location of rotation for dihedral 


F-4 TWIST DEFINITION CARD (7F10-0KRead only If ITW1ST 7D) 



ROUTINE 


GOLUtd 

CODE 

FORMAT 

EXPLANATION 

1 

TWIST Cl) 

Foil 

F10.0 

Twist angle - Dag. for tip airfoil 

11 

TWIST C2) 

Foil 

Twist angle - deg at 2nd airfoil 


F10.0 

station 

21 

TWISTC3) 


Rotation is about L.Eof airfoil. 
Clockwise is positive (l.e trailing 

etc 



edge down) . 


Input 7 twist values per card NS values 
total. Last value Is twist at root 
section.. 



F-5 X DEFINITION CARDS (7F1Q.0) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

XOC(I)* 

Foil 

F10.0 

X/C values used to define airfoil 

11 

etc. 



F-6 INBOARD AIRFOIL 

DEFINITION 

CARD (E10.1)f Input only if NFOILI 0> 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

TYPE I 

Foil 

F10.1 

NACA Airfoil definition either 4 or 
5 digit series, i.e. 0012, 2412, 
64006, etc. 

F-7 INBOARD AIRFOIL UPPER SURFACE DEFINITION CARD fTFlO.OUIttDut only If NF0ILI>0) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

ZIUCI)* 

Foil 

F10.0 

Z/C value for root upper surface 

F-8 INBOARD AIRFOIL LOWER SURFACE DEFINITION CARD (7F10.0) (Input only if NFOILI>0) 


and I VERT = 

0) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

ZILCI)* 

Foil 

F10.0 

Z/C value for root trailing surface 
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F-9 OUTBOARD AIRFOIL DEFINITION CARD (F1Q.1H Input only if NFOILO » 0) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

typed 

Foil 

F10.1 

Same explanation as for TYPE I 


F-10 OUTBOARD AIRFOIL UPPER SURFACE DEFINITION CARD (F10.1) (Input only If NFOILO >0) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

11 

ZOUCD* 

Etc. 

Foil 

F10.0 

Z/C values for tip upper surface 

F-ll OUTBOARD AIRFOIL. LOWER SURFACE. DEFINITION (F10.1) (Input only if NFOILO >0) 


and IVERT 

- 0) 


COLUMN 

core 

ROUTINE 

FORMAT 

EXPLANATION 

1 

ZOL(I)* 

Foil 

F10.0 

Z/C values for tip lower surface 

F-12 SPANWISE PANEL DEFINITION 

CARD (F10.1) (Input only if IND t 0) 

COLUtCI 

core 

ROUTINE 

FORMAT 

EXPLANATION 

1 

ETA(2) 

Foil 

F1Q.0 

Location of span wise sections to be 
defined ETA = y/b/2 


*Input 7 values per card 
NPTS total values 
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SUBROUTINE CONICS 


CO-1 SOURCE PANEL CONTROL CARD (215) (Used on first case only) 

ROUTINE 


COLUMN 

CODE 

FORMAT 


EXPLANATION 

1 

NTOT 

CONICG 

15 


Number of source panels per N-line 

6 

ITRULE 

CONICG 

15 


■ 0 HESS format on output and tape 

» 1 Writes X(I>, Y(I) and Z(I) on 

output and tape 
also =* 0 input more cases 

ft 0 stop 

ITRULE input for 1st case only 

CO-1' is 

the same 

as CO-1 except 

ITRULE 

is not Included. 

CO-2 TRANSLATION 

CARDS (3F10.0, 

2A10) 

(Used on subsequent cases if NTOT <_ 0 return) 

COLUMN 

com: 

ROUTINE 

FORMAT 


EXPLANATION 

1 

XOR 

CONICS 

F10.0 


X Value used to translate all coordinates 

11 

YOR 

CONICS 

F10.0 


Y Value used to translate all coordinates 

21 

ZOR 

CONICS 
FI 0-0 


Z Value used to translate all coordinates 

31 

SECID 

CONICS 

F10.0 


Section ID see HESS input for definition 

41 

XTRPRN 

CONICS 

F10.0 


Print flag used in HESS see HESS 
for details 

CO-3 INITIAL POINT CARD (3F10.0 2A10) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

- 

EXPLANATION 

1 

XE 

CONICS 

F10.0 


X coordinate for 1st point X , in figure 
Co-1 ' J 

11 

YE 

CONICS 

F10.0 


Y coordinate for 1st point Y Q in figure Co-1 

21 

ZE 

CONICS 

F10.0 


Z coordinate for 1st point Z Q in flpure Co-1 



CO-4 CURVE DEFINITION CARD (Al. A9, 15, 5X. FIO.O. 3A1Q) 




ROUTINE 


COLUtti 

com: 

FORMAT 

EXPLANATION 

1 

NMl 

CONICS 

NMI = S Straight curve generated 



A1 

■ C Conic generated 

2 

NM2 

CONICS 

Identifier used so that words straight 



A9 

and curve can be input rather than S and 
C used only for output. 

11 

N 

CONICS 

Number of equally spaced intervals on 



15 

curve generated. 

21 

RHO 

CONICS 

Nondimensioual distance of mid point from 



F10.0 

center of Conic, see figure Co-1 (input 
only if NMl = C) . 

31 

XC 

CONICS 

These are center of Conic, input only if 

41 

YC 

3A10 

Conic does not lie in a primary plant or 

51 

ZC 


Conic is concave otherwise program will 
calculate. 

CO-5 

FINAL POINT CARD (3F10.5. 15) 




ROUTINE 


COLUMN CODE 

FORMAT 

EXPLANATION 

1 

XE 

CONICG 

x coordinate for last point on conic 



F10.5 

(x in figure Co-1) 
n 

11 

YE 

CONICG 

y coordinate for last point on conic 



F10.5 

(y n in figure Co-1) 

21 

ZE 

CONICG 

z coordinate for last point on conic 



F10.5 

(z n in figure Co-1) 

31 

NPARTS 

CONIGG 

number of parts conic is to be divided 



15 

into. This value used only if it is 
desired to vary the point spacing over 
the conic or to output points over more 
or less than the full quadrant. 

CO-6 

CONIC DIVISION CONTROL CARD (215) 

(INPUT ONLY IF NPARTS>0) 



ROUTINE 


COLUMN CODE 

FORMAT 

EXPLANATION 

1 

NJ 

CONICG 

The number of segments on part of the 



15 

conic. NJ is input NPARTS times. The 
order of input is from 1st point to 
last (if data on only part of conic 
desired set NJ*0 on other part not 
desired) . 
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CO-7 CONIC DIVISION VALUES (3F1Q.4. 15) 


COLUMN 

CODE 

ROUTINE 

FORMAT 


EXPLANATION 


1 

XJ 

CONICG 

F10.4 

x value of 

END of segment of 

CONIC 

11 

YJ 

CONICG 

F10.4 

y value of 

end of segment of 

conic 

21 

ZJ 

CONICG 

F10.4 

z value of 

end of segmznt of 

conic 


NPARTS number of these values input XJ, YJ, ZJ must progress 
on along conic line from first to last point (for part of 
conic not desired set starting and end points the same 
except for slight drift to prevent numerical bomboffs). 


*o» Xo. z o 



RHO 


Dm 

Dl 


Figure CO-1 
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SUBROUTINE CIRCGN 


This subroutine generates data for inlets - half inlets - circle - half circle. 


C-l PROGRAM CONTROL CARD (2A4. A2. 3X. 12. 5X. 2A10) 


COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

N1 

CIRCGN 

A4 

N1 * half do half circle or half inlet 

N1 * circ do circle 

N1 * inle do inlet 

= F3D stops execution CTTC'^T 

5 

N2 

CIRCGN 

A4 

used to complete definition of surface to 
consider 

if N1 * half then 

N2 — INL do half inlet 

N2 *-CIR do half circle 

o 

N3 

CIRCGN 

A2 

used to complete definition of surface, 
used in output only 

14 

N 

CIRCGN 

12 

input for circles and half circles only. 

The number of divisions of the angular 
location of pts (&$> 2v (circles) 

N 

half circles) 

N 

Section ID - Hess in^ut - see hess for details 

21 

LI 

CIRCGN 

A10 

il 

Note: 

L2 CIRCGN Print Flag - See HESS input for details 

A10 

»1 + N2 + N3 - inlet 

half-inlet 

circle 

half-circle 

These are only acceptable inputs to N1 N2 & NG 

C2 CURVE DEFINITION CARD (4F10.0. 

A10) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

0 

CIRCGN 

F10.0 

Diameter of circle on inlet (input only if R*0) 

11 

R 

CIRCGN 

F10.0 

Radius of circle or inlet (input only if D*0) 

21 

PHIZ 

CIRCGN 

F10.0 

Initial value of angular position input for 
circles or half circles only 

31 

PH1F 

CIRCGN 

F10.0 

Final value f angular position input for circles 
or half uire us only 



C-3 CURVE 

ROTATION CARD (4F10.0. A10) 


CODE 

ROUTINE 

FORMAT EXPLANATION 

1 

XC 

CIRCGN x center of rotation 

F10.0 

11 

YC 

CIRCGN y center of rotation 

F10.0 

21 

ZC 

CIRCGN z center of rotation 

F10.0 


31 

ROT 

CIRCGN 

F10.0 

rotation about pitch axis about (x c , y c , z fi ) 
degrees 

41 

N0ENDFL 

CIRCGN 

F10.0 

if blank ENDFILE written on tape 
if nonblank no ENDFILE written 

when done with CIRCGN input N1 * end 
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SUBROUTINE HESPLT 


Subroutine BESPIT generates CALCOMP plots of the bess input geoaetry. The 
HESPLT input consists of all of the hess input plus sane viewing angle input 
data. Therefore we will only show the extra data required here and will return 
the user to the hess input section for the nain body of Input data. 

HESS DATA. INPUT (See hess input) 

Following HESS input the following card or cards are added 


HP-1 VIEWING AHGLE COBTROL CARDS (4F10.3.20X.2A10) 




ROUTINE 

■ 

maw 

CODE 

FORMAT 

EXPLANATION 

1 

ALPHA 

HESPLT 

pitch angle at desired view deg (+ left hand rule) 
if ALPHA input greater than 360° then program 
terminates . Use this option to END DATA input. 

li 

BETA 

HESPLT 

F10.3 

PAW angle of desired view (POS according to right 
hand rule) 

21 

PHI 

HESPLT 

F10.3 

Roll angle at desired view (positive according 
to left hand rule) 

31 

<nrar.B 

HESPLT 
no. 3 

Plot scale 

If scale » 0 the routine will determine the proper 
scale. 

61 

VW(1) 

HESPLT 

A10 

to character alpha numeric array describing view 

71 

VW(2) 

HESPLT 

A10 

(i.e. top view) 


Input as nany of the viewing angles as desired at this point 
stopping input with a value >360 in the alpha input array. 
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SUBROUTINE JUNC PAM 


This routine finds Che juncture between a nonlifting HESS section (fuselage, 
nacelle) and a lifting section (wing, horizontal, vertical). The input 
consists of coordinate data set up for hess input. The data required is listed 


below. 

For details 

on the proper order of data see the HESS input module. 

TITLE CABO (SAlO) 



COLOMf 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

TELE (I) 

JUNC 

8A10 

Title of case 

HON LIFTING SECTION CONTROL CARD (1613) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

HLSTR 

JUNC 

15 

Nuaber of non-lifting strips in non-lifting section 

6 

NLS0R 

JUNC 

15 

Number of panels per strip in non lifting section 

LIFTING SECTION CONTROL CARD (1615) 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

NSTRIP 

JUNC 

15 

Nuaber of lifting strips in lifting section 

6 

NS0RCE 

JUNC 

15 

Number of panels per strip in lifting section 

11 

NWAKE 

JUNC 

15 

Number of wake panels per strip in lifting section 

NON LIFTING SECTION 

INPUT CARDS 

(T61.A10.Tl t (6F10.4)} 

COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

XNL(I) 

JUNC 

F10.4 

x coordinates on 1st N line of non lifting section 

11 

YNL(l) 

JUNC 

F10.4 

y coordinate on 1st N line of non-lifting section 

21 

ZNL(I) 

ZNL(I) 

JUNC 

F10.4 

z coordinate on 1st N line of non- lifting section 
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COLUMN 

ROUTINE 

CODE FORMAT 


31 

XNL(I+) JCNC 
F10.4 


41 

YNL(I+1) JUNC 
F10.4 


51 

ZNL (1+1) JUNC 
F10.4 


61 

SECIDNL JUNC 
A10 

Section Identifier 


z sets of ZNL, 

TNL, ZNL per card Input NLS0R+1 values 


{6F10.4} 

COLUMN 

ROUTINE 
CODE FORMAT 

EXPLANATION 


1 

XNL(I) 

JUNC 

F10.4 

x coordinate on N line of non lifting coordinates 

11 

TNL (I) 

JUNC 

F10.4 

y corresponding to x above 

21 

ZNL (I) 

JUNK 

F10.4 

z corresponding to x above 

31 

XNL(I) 

JUNK 

F10.4 



etc 


Input 2 sets of XML, YML, ZNL per card. 

Input NLSDR+1 values. 

Repeat these cards for all M lines In this section (1*2, NLSTR+1) . 


LIFTING SECTION INPUT CARDS (T61.A10.T1. (6F10.4)) 




ROUTINE 



COLUMN 

CODE 

FORMAT 

EXPLANATION 

1 

XL(I) 

JUNC 

F10.4 

x coordinate on 

1st N-llae of lifting section. 

11 

TL(I) 

JUNC 

F10.4 

y coordinate of 

1st N-line of lifting section 

21 

ZL(I) 

JUNC 

F10.4 

z coordinate of 

1st N-line of lifting section 
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COLUMN 


CODE 


ROUTINE 

FORMAT 


EXPLANATION 


31 

XL(I+) 

JCNC 

F10.4 



41 

TL(I+1) 

JUNC 

F10.4 



51 

ZL(I+1) 




61 

SECIDL 

JUNC 

A10 

Lifting Section Identical 



Input 2 sets of XL, YL, ZL per 
Input NSORCE + NWAKE+1 values 

card 

WML 

.♦i in * o: mm ah ' a * mai 



COLUMN 

CODE 

ROUTINE 

FORMAT 

EXPLANATION 

1 

XL(I) 

JUNC 

F10.4 

z coordinate 

on Nline of lifting section 

11 

YL(I) 

JUNC 

F10.4 

y coordinate 

on N line of lifting section 

21 

2L(I) 

JUNC 

F10.4 

z coordinate 

on N line of lifting section 

31 

XL(I) 

JUNC 

F10.4 

x coordinate 

on N line of lifting section 


Input 2 sets of XL, YL, ZL per card. 

Input NSORCE + NWAKE -t-1 values for each N line. 
Repeat input NSTRIP + 1 times. 
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APPENDIX 3 - INTERACTIVE GRAPHICS INPUT PROGRAM (HESTER) 


The input geometry required by the Hess three-dimensional program 
must be developed carefully in order to adequately describe the con- 
figuration being analyzed. Most aircraft configurations are complex 
and thus require a considerable manhour expenditure to develop an ap- 
propriate model. For example, the modeling of a typical V/STOL configura- 
tion, without using any of the input routines discussed in this report, 
could take two to four weeks of concentrated effort. Even using the 
routines in this report, considerable effort is still required and the 
possibility of mistakes in the geometry is high due to the large number 
of coordinates involved. Therefore in order to reduce both the model 
development time and the possibility of errors, a stand alone interactive 
graphics program was developed. This program permits the user to display 
all or part of the model developed and to make corrections where necessary. 
The program has the capability of moving points, lines or sections, adding 
or deleting lines and deleting sections. Note that individual points may 
only be moved, they can not be added or deleted. 

The following discussion describes the operation of this program. 
PROGRAM DISPLAY CAPABILITIES 


When program HESTER is accessed, a graphical representation. Figure 
B-l, will appear on the screen. The current values of scale factor and 
pitch, roll and yaw angles will be displayed. This graphical display can 
be altered in several ways as explained below. 

A) PLOT SCALING 

The size of the graph being displayed can be changed in several 

ways. 

(1) The size can be increased by a factor of 2 by pressing the 
(>) (greater than) key. The display will erase and be re- 
plotted twice the original size. This process can be re- 
peated as many times as desired. 

(2) The size can be decreased by a factor of 2 by pressing the 
(<) (less than) key. The display will erase and be replotted 
at one-half the original size. This process can be repeated. 

(3) The scale parameter can be changed to whatever the user 
desires by pressing the (I) key. The program will then 
respond with a request and the scale desired can be input. 

B) PLOT POSITION 


In the process of scaling, it sometimes becomes necessary to re- 
position the viewing window so that a particular region of the 
configuration is at the center. This is accomplished by moving 
the cross hairs to the point on the current graph that is desired 



to be in the center. The (*) (equal) sign key is now depressed 
which will cause Che graph co be replotted at the same scale with 
this new position in the center of the graph. 

C) CONTENTS OF DISPLAY 


The contents of the display may also be altered. It may be desired 
to only view one section or possibly a group of sections alone. 

This can be done by hitting the (S) key. The program will respond 
by asking what section or sections are desired. The sections are 
selected either by section number, section ID, or section type, 
(i.e., all non-lifting sections may be plotted) (Figure B-2). 

PROGRAM MODIFICATION CAPABILITIES 


Once a graphical representation of the model is available, it sometimes 
becomes desirable to modify the model by: (1) moving points, lines, or even 

sections, (2) adding or deleting a strip-line or panel-line, and (3) delet- 
ing a section.. The programs contain several options as described below for 
performing these tasks. 

A) MOVING A POINT. LINE OR SECTION 

Either an individual point, a particular line or an entire section 
may be moved as explained below. 

(1) MOVEMENT OF A POINT 

Position the crosshairs on the point to be moved and depress 
the (.) (period) key. A slash will then be placed on this 
point for identification. Then move the cursur crosshairs 
to the position desired and depress the (@) key. The plot 
will be redrawn with the point moved. 

(2) MOVEMENT OF A LINE 

If an entire line is to be moved, a similar procedure is 
followed. Position the cursur crosshairs over one point 
in the line, at an intersection point, and depress the (L) 
key. Then, in order to identify the line, move to another 
point on the line and again depress the (L) key. Note, if 
desired several points on the line may be chosen as above. 
Then, position the cursur crosshairs where it is desired to 
move the first point on the line, and depress the (@) key. 

The entire line will then be moved to this new location. 

This procedure works on either a strip-line or a panel-line. 

(3) MOVEMENT OF A SECTION 

An entire section may be repositioned if desired, by placing 
the cursur crosshairs at a control point selected within the 
section and depressing the (T) key. This control point must 
be a point for which the transformed location is known. Then 
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move the cursors to the location where the control point is 
to be transformed and depress the ((?) key. The plot will be 
redrawn with the section containing the control point moved 
to the new location. 

B) ADDING A LINE 

A line may be added to a configuration by the following procedure: 

(1) position the cursur at a point where a line is desired and de- 
press the (.) (period) or (L) key. Then move the cursur to another 
point where the line is pass through and depress the (+) key. A 
line will then be formed from the existing data by interpolation 
which passes through the two points chosen above. 

C) DELETING A LINE OR SECTION 


A line or. section may be' deleted as explained below 

(1) DELETING A LINE 

Place the cursur on one point in the line and press the (.) 
(period) or (L) key. Then move the cursur to another point 
on the line and depress the (-) key. This will remove that 
line from the section. 

(2) DELETING A SECTION 

Place the cursur on a point in a section and depress the (T) 
key. Then depress (-) key to delete section. 

PROGRAM EDITING CAPABILITIES 


The data set from which the plots are drawn may be edited as explained 
below. The data set can be viewed by depressing the (E) key. This causes 
the data set to be printed on the screen one page at a time. The following 
procedure is used to modify values in the data set. 

(1) CHANGING AN ENTRY 

To change a value in the data set, move the cursur to the 
value in question and depress the (C) key and then enter 
the new value. This is repeated as needed. 

(2) VERIFICATION OF CHANGES 

After making changes, depress the (V) key, this will cause 
the same page to be rewritten with the changes incorporated. 

(3) LISTING OF DATA 

After a change in made or after a page has been reviewed, 
depress the (g) key to continue listing data. 
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(4) 


SKIPPING DATA 


If it is desired to skip parts of a section, simply de- 
press the (N) key. This will cause the program to skip 
to the next data section. 

(5) EXITING FROM EDIT MODE 

Depressing the (:) (colon) key causes the program to return 
to the program mode. 

There are several other operations which may be performed when operating 
this program. These operations are discussed below. In addition, the loca- 
tion of the modified data files are discussed. 

(1) At any time in the modification of a data set, the current 
data may be saved by depressing the (C) key. .This creates 
a file called CHKDAT. Everytimie the (C) key is pressed, 
the file is rewritten. 

(2) 7i mistakes are made in a change and it is desired to return 

to the previous data set, simply depress the (N) key. This 

deletes the current change. 

(3) If it is desired to return to the original data set depress 
the (0) key. 

(4) If it is desired to view the file saved or the CHKDAT file, 
depress the (R) key. 

Finally, when all modifications and operations are finished, depress the (B) 
key. This terminates the program. At termination, three files are available: 

(1) the original data file is available under the input file name and tape 5; 

(2) the CHKDAT file is available and (3) the final file which was being worked 

on at termination is stored on the NEWDAT file. Note tht by proper use of the 

CHKDAT file, it is possible to create two files which contain different modi- 

fications to the same configuration. 
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ORIGINAL PAGE 19 
OF POOR QUALITY 


SCALE « 1.0 
PITCH = 0.0 
YAW = 0.0 
BOLL = 90.0 



Figure B-l. Original Plot placed on screen 
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ORIGINAL PAGE H 
OF POOR QUALITY 


SCALE = 1.0 

PITCH = 0.0 

YAW * 0.0 

POLL = 90.0 

ALL 

NONLIFTING 

LIFTING 

INLETS 

1 MAIN SECT 

? SECTION 2 

3 SECTION 3 

4 SECTION 4 



Figure 3-2. Plot displayed when (S) Key Depressed 
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